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In this review, we discuss the confining and finite-temperature properties of the 3D 

SU(Af) Georgi-Glashow model, and of 4D compact QED. At zero temperature, we 

O , ' derive string representations of both theories, thus constructing the SU(Af)-version of 

^ , Polyakov's theory of confining strings. We discuss the geometric properties of confining 

,.S^ i strings, as well as the appearance of the string 6-term from the field-theoretical one in 

4D, and fc-string tensions at A'^ > 2. In particular, we point out the relevance of nega- 
tive stiffness for stabilizing confining strings, an efi^ect recently re-discovered in material 
J^ ' science. At finite temperature, we present a derivation of the confining-string free en- 

\-i ergy and show that, at the one-loop level and for a certain class of string models in the 

^ large-D limit, it matches that of QCD at large A'^. This crucial matching is again a con- 

sequence of the negative stiffness. In the discussion of the finite-temperature properties 
of the 3D Georgi-Glashow model, in order to be closer to QCD, we mostly concentrate 
at the effects produced by some extensions of the model by external matter fields, such 
as dynamical fundamental quarks or photinos, in the supersymmetric generalization of 
the model. 
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3D Georgi-Glashow model and confining strings 3 

Many papers related to this review would have never been written without 
Ian Kogan's seminal ideas on the finite-temperature phase transitions in 
string theory and in the 3D Georgi-Glashow model. The sudden death of 
Ian Kogan was a loss for physics and a personal tragedy for everybody who 
had known him, in particular for the authors of this manuscript. In the 
personal sense, Ian was a very bright man with a warm heart, and he will he 
deeply missed as a good friend. We nevertheless strongly believe that Ian 's 
ideas will survive over decades, inspiring many new generations of physicists. 

1. Introduction 

During last 30 years, the problem of quark confinement remains as a great 
challenge not only to QCD theorists, but to the whole theoretical high-energy 
physics community. In general, confinement can be defined as the absence 
in the spectrum of a certain field theory of physical |in) and |out) states 
of some particles, whose fields are nevertheless present in the Lagrangian 
of that theory (see e.g. [1~3] for reviews and books on confinement). With 
regard to QCD, this definition means the absence of asymptotic quark and 
gluon states, i.e. states which carry color. '^ This fact is refiected in the 
linear growth with the distance of the potential between two color parti- 
cles, as well as in the logarithmic growth of the strong-coupling constant. 
The latter makes the standard perturbative diagrammatic techniques inap- 
plicable at distances larger than O ( Aq(^p j. However, such distances are 
of primary interest, since only there physically observable colorless states 
of hadrons are formed, whereas at smaller distances one deals with the un- 
physical colored states - quarks, gluons, and ghosts. A question, which may 
naturally be posed at this point, is whether the QCD Lagrangian, operating 
with these colored states only, eventually yields the correct description of 
colorless degrees of freedom as well. A reliable indication that this is really 
the case comes from the simulations of the QCD Lagrangian on a lattice 
(for recent reviews see e.g. [2,5,6]), which yield a reasonably good descrip- 
tion of hadronic spectra. However, as of today, a systematic analytic way 
of describing large-distance effects in QCD in terms of microscopic (colored) 
degrees of freedom is unfortunately lacking. The breakdown of the perturba- 
tive expansion at the large distances under discussion naturally introduces 
the notion nonperturbative for these effects, as well as for the techniques 



^ An attempt to construct a momentum-space interpretation of tiiis phenomenon lias been done 
(see e.g. Ref. [4]). In this review, we will, however, prefer to deal with the conventional space-time 
picture, which enables one to directly operate with such notions as the vacuum correlation length 
or potential between color particles. 



February 1, 2008 18:47 WSPC/Trim Size: 9.75in x 6. Sin for Proceedings Antonov 



4 D, Antonov and M. C. Diamantini 

which physicists attempt to invent for their description. 

It is actually not a great fault that the perturbative expansion is inappli- 
cable at large distances, since, being formulated in terms of colored states, 
it is anyway not suitable to keep gauge invariance under control. Indeed, 
an individual diagram is always gauge-dependent, merely because a certain 
gauge fixing should be performed before one starts to compute diagrams. 
Contrary to that, hadronic states, being colorless, are gauge-invariant. More- 
over, since the QCD vacuum possesses an unbroken gauge symmetry,^ only 
gauge-invariant amplitudes, being averaged over the vacuum, yield a nonva- 
nishing result.'^ Although the gauge dependence disappears when one sums 
up contributions of individual diagrams to a certain gauge-invariant quan- 
tity, it would be more natural to have an expansion already operating with 
such quantities at any intermediate stage. Such an expansion, which is the 
expansion in the inverse number of colors, has been proposed by 't Hooft [9] 
and further developed in the framework of loop equations [10] (see e.g. [11] 
for reviews). This is a classic example of a nonperturbative approach to 
QCD. 

Another nonperturbative phenomenon in QCD, which is of the same fun- 
damental importance as confinement, and whose complete analytic under- 
standing is also still lacking, is the spontaneous chiral symmetry breaking 
(SCSB). This is the symmetry of the QCD Lagrangian with Nf massless 
flavors under the global transformations, which are the U(iVj)x U(A'^^) in- 
dependent rotations of left- and right-handed quark fields. These are equiv- 
alent to the independent vector and axial U(A'^^) rotations of the full four- 
component Dirac spinors, under which the QCD Lagrangian remains invari- 
ant as well. At the same time, the axial transformations mix states with dif- 
ferent P-parities. Therefore, were the chiral symmetry unbroken, one would 
observe parity degeneracy of all the states whose other quantum numbers 



^ We do not discuss in this review tiie modern line of research devoted to the so-called color 
superconductivity [7]. In the color superconducting phase of QCD, which takes place at high 
baryonic densities corresponding to the values of the quark chemical potential > 0(400 MeV), 
the gauge symmetry is spontaneously broken. 

"^ To get a description of large-distance effects in terms of microscopic degrees of freedom, it looks 
natural to integrate over the high-momentum modes in the same way as is usually done in a 
derivation of renormalization-group (RG) equations in statistical mechanics. For instance, this 
approach gives the correct RG behavior of the running strong-coupling constant (at the one-loop 
level) [3]. However, it typically suffers from the above-mentioned problem of violation of gauge 
invariance, as far as RG equations for Green functions are concerned. Recent progress in this 
direction of research is concentrated around the concept of the so-called exact RG flow (see e.g. 
Ref. [8] for a recent review). 
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3D Georgi-Glashow model and confining strings 5 

are the same. The observed sphttings between such states are, however, too 
large to be explained by the small quark masses. Namely, this splitting is 
of the order of hundreds of MeV, whereas the current masses of light u- 
and d-quarks are of the order of a few MeV.'^ This observation tells us that 
the chiral symmetry of the QCD Lagrangian is broken down spontaneously. 
The phenomenon of SCSB naturally leads to the appearance of light pseu- 
doscalar Goldstone bosons, whose role is played by pious, which are indeed 
the lightest of all the hadrons. The order parameter of SCSB is the chiral 
quark condensate {"ipip) ^ — (250MeV)^. Its appearance is quite natural in 
QCD, where, due to the strong attraction between quarks and antiquarks, 
it costs very little energy to create a massless quark-antiquark pair. Having 
zero total momentum and angular momentum, such pairs carry net chiral 
charge, hence is the notion "chiral condensate". 

Since the early days of QCD, when it has been realized that its nonpertur- 
bative phenomena, confinement and SCSB, are of fundamental importance, 
numerous nonperturbative approaches have been proposed in an attempt to 
describe one or other of these two phenomena in a controllable way.*^ On 
the side of confinement, these approaches include e.g. the already men- 
tioned expansion in the large number of colors [9], loop equations [10], the 
stochastic vacuum model [1], method of Abelian projections [14]. On the 
side of SCSB, the classic approach is the one based on the Nambu-Jona- 
Lasinio models [17], as well as on the related nonlinear chiral meson La- 
grangians [18] (the latter approach has further been developed in [19]; for a 
review see [20]). There also exist microscopic models of SCSB in QCD, based 
on instantons [21] or dyons [22], which have been put forward in Refs. [23] 
and [24], respectively. 

Although a lot of physical effects have been captured by the above- 
mentioned methods, unfortunately none of them provides the full solution 
of QCD. The concept of "full solution" has several facets. According to the 



'^ The current mass of the s-quark, which is around 150 MeV, is still smaller than the typical 
splitting values by at least a factor of three. 

"^ Unfortunately, practically no approach exist, which describes both phenomena. Partially, that is 
because it is not yet fully accepted that these phenomena are related to each other microscopically, 
i.e. that the same vacuum configurations are responsible for both phenomena. For instance, as 
has been shown in Ref. [12], this is not so in the case of the so-called Faddeev-Niemi effective 
action [13]. An example of a nonperturbative approach, where the two phenomena are related 
to each other, is the stochastic vacuum model of QCD [1], which does not refer to a particular 
microscopic vacuum configuration. 

' A comprehensive lattice analysis of the latter has recently been performed [15] (see Ref. [16] for 
reviews). 
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conventional understanding, the solution of a certain field theory means a 
prescription of how to compute an arbitrary gauge-invariant vacuum ampli- 
tude. In QCD, the results should be in an agreement with the experimental 
and lattice data on the properties of hadrons. The gluonic and quark aver- 
ages should then give a correct quantitative description of confinement and 
SCSB, respectively. The values of local averages (condensates) and relations 
between these quantities should also agree with those known from QCD 
sum rules [25]. Furthermore, the standard diagrammatic expansion should 
be reproducible (giving, in particular, asymptotic freedom), and the non- 
perturbative dimensionful quantities should be made expressible in terms of 
the only dimensionful QCD parameter, Aqcd- Moreover, some constraints, 
which are accepted to be rigorous in QCD, should be obeyed. For instance, 
the above-mentioned expansion in the large number of colors should be cor- 
rectly reproducible and should respect the large- A^ loop equation. The so- 
lution of QCD should also accommodate classical vacuum configurations of 
the gluodynamics action, such as instantons. 

In this review, we are going to discuss the 3D Georgi-Glashow (GG) 
model, which is the QCD-related model possessing the property of confine- 
ment [26] . Our primary goal will therefore be the study of confinement (and 
not of the SCSB), as well as of the deconfining phase transition at finite 
temperatures. The advantage of the 3D GG model with respect to QCD is 
that confinement in it takes place in the weak-coupling regime. It turns out 
that, already in this regime, the vacuum (i.e. the ground state) of this model 
is nonperturbative, being populated by 't Hooft-Polyakov monopoles [27], 
which provide the permanent confinement of external fundamental charges. 
As a guiding principle of our analysis we will use the string picture of con- 
finement, therefore let us briefly discuss it. 

In QCD, the linearly rising confining interquark potential is associated 
to a string-like configuration of the gluonic field between quarks, usually 
called the QCD string. Indeed, the energy of a string grows linearly with 
its length, E{R) = aR .^ According to the Regge phenomenology, the string 
energy density o", called the string tension, is approximately (440 MeV)^. 
The string can naturally be called confining (the notion, which is always 
used in confining gauge theories other than QCD, such as the 3D GG model), 
since with the increase of the distance R, it stretches and prevents a quark 



s In this review, we only briefly discuss the phenomenon of string breaking (see the end of sub- 
section 4.2). String breaking always happens at a certain distance if dynamical matter fields, 
transforming by the same representation of the gauge group as the confined external ones, are 
present. 
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3D Georgi-Glashow model and confining strings 7 

and an antiquark from the separation to macroscopic distances. As for any 
dimensionful quantity in QCD, the string tension should be proportional to 
the respective power of Aqcd • Namely, 



o" « Aqcd = a ^ exp 






a exp 



167r2 



{^^N,-lNf)g^{a-^) 



XI) 



where a ^ stands for the inverse UV cutoff (e.g. the lattice spacing) 
Furthermore, "~" means "at the one-loop level", at which 



*^'--(y"-Hil^ 



As is seen explicitly from Eq. (1.1), all the coefficients in the expansion of 
a in (positive) powers of g^ vanish, which means that the QCD string is 
indeed an essentially nonperturbative object. The string should therefore 
be produced by some nonperturbative background fields. On top of these, 
however, one expects to have some quantum fluctuations of the gauge field, 
which give rise to the string excitations. 

The confining quark-antiquark potential corresponds to the so-called area 
law of the Wilson loop: ^ 



{W{C)) ^ — ( trP exp ( ^g j> A^T^dx, | ) '"^ e 



Here, Sjnin(C) is the surface of the minimal area, bounded by the trajectory 
C of the quark-antiquark pair, and | . . . | means either a length or an area. 
The confining string, which sweeps out the minimal surface, can naturally 
be viewed as a product of the above-mentioned strong background fields. In- 
stead, quantum fluctuations around these enable the string to sweep out with 
a nonvanishing probability any other surface S(C), different from Smin((7). 
To derive the string representation of QCD would be to give sense to the 



It should be compared with the perimeter law, 



which corresponds to the Coulomb potential and is found in non-confining theories, e.g. (non- 
compact) QED. 
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formula 

(VF(C)) ^ (H^(S„,i„(C))) = Y, e-^P(^)l . (1.2) 

S{C) 

Here, Yl ^^^ ^i^i^)] stand for a certain sum over string world sheets and 

S(C) 
a string effective action, both of which are yet unknown in QCD. 

The first problem one encounters in trying to determine S'[S(C)] is that 
fundamental strings [28] can be quantized only in critical dimensions: a con- 
sistent quantum theory describing strings out of the critical dimensions has 
not yet been found. The simplest model, the Nambu-Goto string, can be 
quantized only in space-time dimension D = 26 or D < 1 because of the con- 
formal anomaly. It is appropriate to describe an effective string picture for 
confinement in QCD , but is inappropriate to describe fundamental smooth 
strings dual to QCD [29], since large Euclidean world sheets are crumpled. 
The picture of a fundamental string theory dual to QCD is strongly sup- 
ported by a recent lattice calculation by Liischer and Weisz [30], where evi- 
dence of a string behavior in the static quark-antiquark potential has been 
found down to distances of 0.5 fm. 

In the rigid-string action [31,32], the marginal term proportional to the 
square of the extrinsic curvature, introduced to avoid crumpling, turns out 
to be infrared irrelevant and, thus, unable to provide smooth surfaces. 

Recent progress in this field is based on a new type of action. In its 
local formulation [33,34], the string action is induced by an antisymmetric 
tensor field. This action realizes explicitly the necessary zig-zag invariance 
of confining strings [34,35]. It can be derived without extra assumptions [36] 
for the confining phase of compact U(l) gauge theories [3]. An alternative 
approach to the induced string action was originally proposed in [37], and is 
based on a five-dimensional, curved space-time string action with the quarks 
living on a four-dimensional horizon [38]. The formulation of the string 
theory in the five-dimensional curved space-time is closely related to the 
AdS/CFT (Anti de Sitter/Conformal Field Theory) correspondence [39]. In 
fact, with a special choice of the metric in the curved space, one recovers 
the AdSs space, thereby providing a string theory description of a conformal 
gauge theory [39]. 

The main characteristic of the effective string action obtained by inte- 
grating out the tensor field is a non-local interaction with negative stiffness, 
that can be expressed as a derivative expansion of the interaction between 
surface elements. To perform an analytic analysis of the geometric proper- 
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ties of these strings, this expansion must be truncated: this clearly makes 
the model non-unitary, but in a spurious way. Moreover, since the stiffness 
is negative, a stable truncation must, at least, include a sixth-order term in 
the derivatives. The role of negative stiffness, as first pointed out in [36,40] 
is crucial. It is in fact the sixth-order term, forced by the negative stiffness, 
that suppresses the formation of spikes on the surfaces and leads to a smooth 
surface in the large-D approximation. In fact, in [40,41] it has been shown 
that, in the large-D approximation, this model has an infrared fixed point 
at zero stiffness, corresponding to a tensionless smooth string whose world 
sheet has Hausdorff dimension 2, exactly the desired properties to describe 
QCD flux tubes. As first noticed in [42,43], the long-range orientational or- 
der in this model is due to an antiferromagnetic interaction between normals 
to the surface, a mechanism confirmed by numerical simulations [44]. The 
presence of the infrared fixed point does not depend on the truncation [41] 
and it is present for all ghost- and tachyon- free truncations. Moreover, the 
effective theory describing the infrared behavior of the confining string is a 
conformal field theory with central charge c = 1. 

Another important feature of the negative-stiffness model is its high- 
temperature behavior. Contrary to all previous string models for QCD, it is 
able to reproduce the large- A^ QCD behavior, found by Polchinski and Yang 
in [89], in both sign and reality properties [45]. 

It is remarkable that the role of negative stiffness, while first discovered 
in the context of strings [36,40] and membranes [47], has been rediscovered, 
and actually experimentally tested, also in material science [48]. In fact, 
it has been found that composites with negative stiffness inclusions have 
higher overall stiffness than that of their constituents. Such composites find 
applications in which high stiffness and damping are needed, permitting 
extreme properties not previously anticipated. 

Equation (1.2) would clearly be a generalization of a path-integral rep- 
resentation for a propagator of a particle, for example boson, (0(x)0(x')) = 
^ e~^^^xx'\ where P^x' is a path connecting the points x and x'. In par- 

ticular, within this analogy, the role of the classical trajectory of a particle 
would be played by Sjnin(C). For a point-like particle, the measure in the 
sum over paths is known and depends only on the dimensionality of the 
space-time. The world-line action ^[P^:!:'] can also be evaluated, either ana- 
lytically [for certain potentials V{(p) or external gauge fields if (p is charged], 
or using Feynman's variational method. On the string side, a derivation of 
the measure in the sum over world sheets has been discussed in Ref. [49] 
for the case of the Abelian Higgs model in the London limit. As will be 
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seen below, in a certain case the string effective action derivable in the 3D 
GG model is the same as that of the London limit of the dual Abelian Higgs 
model. The summation over string world sheets is, however, realized in these 
two models in different ways. Namely, in the dual Abelian Higgs model it 
stems from the integration over singularities of the dual Higgs field (vortex 
cores), which is already present in the original partition function. Instead, 
the 3D GG model does not contain any dual Abelian Higgs field, and the 
summation over world sheets in this model is realized by means of the result- 
ing string effective action itself [34] . Numerous investigations of this action 
have further been performed (see e.g. Refs. [36,42,43]), and some of these 
will be discussed in this review. A separate section will be devoted to the 
finite-temperature properties of confining strings [45]. 

For the purpose of the study of the so-called fc-strings and merely to be 
closer to real QCD, we will deal with the SU(A^)-generalization of the stan- 
dard SU(2) GG model, which will be introduced in the next section. As for 
the (A'^ = 2) 3D GG model, it is a classic example [26] of a theory which 
allows for an analytic description of confinement. As has already been men- 
tioned, confinement in this model is due to the plasma of point-like magnetic 
monopoles, which produce random magnetic fluxes through the contour of 
the Wilson loop. In the weak-coupling regime of the model, this plasma is 
dilute, and the interaction between monopoles is Coulombic, being induced 
by the dual-photon exchanges. Since the energy of a single monopole is a 
quadratic function of its flux, it is energetically favorable for the vacuum 
to support a configuration of two monopoles of unit charge (in the units of 
the magnetic coupling constant, g-m), rather than a single, doubly-charged 
monopole. Owing to this fact, monopoles of unit charge dominate in the 
vacuum, whereas monopoles of higher charges tend to dissociate to them. 
Summing over the grand canonical ensemble of monopoles of unit charge, 
interacting with each other by the Coulomb law, one arrives at an effective 
low-energy theory, which is a 3D sine-Gordon theory of a dual photon. The 
latter acquires a mass (visible upon the expansion of the cosine potential) by 
means of the Debye screening in the Coulomb plasma. The appearance of 
this (exponentially small) mass and, hence, of a finite (albeit exponentially 
large) vacuum correlation length is crucial for the generation of the funda- 
mental string tension, i.e. for the confinement of an external fundamental 
matter. It is worth noting that a physically important interpretation of these 
ideas in terms of spontaneous breaking of magnetic Z2 symmetry has been 
presented in reviews [50] and Refs. therein. 

While the confining properties of the 3D GG model have been well 
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known since Polyakov's pioneering paper [26], the finite-temperature prop- 
erties of this model have been addressed only recently, starting with the 
papers [51,52]. It turns out [52] that charged matter fields of W-bosons 
play the crucial role for the dynamics of the phase transition. Below we will 
review this issue and also discuss the influence of other matter fields on the 
finite-temperature properties of the model. Such fields are either already 
present in the original Lagrangian (e.g. Higgs [53], or dual photinos in the 
supersymmetric generalization of the model [54]), or can be included in the 
framework of a certain extension of the model (e.g. massless fundamental 
quarks [55] or heavy fundamental bosonic matter [56]). 

The outline of this review is the following. In the next section, we will 
introduce the 3D GG model in the general SU(A'^) case. In section 3, we 
will find the string tension of the fundamental Wilson loop defined at a flat 
contour (henceforce called for shortness "flat Wilson loop"). In section 4, we 
will develop a theory of conflning strings based on the Kalb-Ramond field, 
which enables one to deal with non-fiat Wilson loops. Using for concreteness 
the case of fundamental representation, we will present in subsection 4.1 two 
methods by means of which the theory of confining strings can be derived. In 
subsection 4.2, the case of the adjoint Wilson loop will be considered, and the 
corresponding theory of confining strings will be constructed in the large- A^ 
limit. In subsection 4.3, we will study the spectrum of /c-strings, i.e. strings 
between sources in (higher) representations with a nonvanishing A^-ality. 
These sources carry a charge k with respect to the center of the gauge group 
Z]\f and can be seen as a superposition of k fundamental charges. Clearly, 
the spectrum of such strings is an important ingredient for the complete 
description of the confining dynamics of the 3D GG model. In section 5, the 
SU(A^)-theory of confining strings will be generalized to the SU(iV)-version 
of 4D compact QED (in the continuum limit) with the field-theoretical 9- 
term. As has been found in Refs. [33,43], for the usual compact QED, 
this term leads to the appearance of the string 0-term. The latter, being 
proportional to the number of self-intersections of the world sheet, might 
help in the solution of the problem of crumpling of large world sheets [3,31]. 
The critical values of 6, at which this happens, will be derived in the general 
SU(A'^)-case for fundamental and adjoint representations, as well as for k- 
strings. In section 6, various geometric features of confining strings will 
be studied. In section 7, the thermodynamics of confining strings will be 
discussed, and a derivation of the one-loop free energy of a string in the 
large-D limit will be presented. Again, we will show that it is the presence 
of negative stiffness that allows one to reproduce the large-A'^ behavior of 
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high-temperature QCD. 

In section 8, we wih pass from the thermodynamics of confining strings to 
the thermodynamics of the 3D GG model itself. After an introduction to this 
subject in subsection 8.1, we will pay particular attention to the influence of 
various matter fields on the dynamics of the deconfining phase transition. In 
subsections 8.2 and 8.3, we will consider an approximation where W-bosons 
are disregarded. Subsection 8.2 will be devoted to the influence of the Higgs 
field (when it is not infinitely heavy) on the RG flow, while in subsection 8.3 
we will consider the model in the presence of external dynamical fundamental 
quarks. In subsection 8.4, we will first discuss the crucial role of W-bosons 
in the dynamics of the phase transition in the finite-temperature 3D GG 
model and then consider the supersymmetric generalization of the model. 
In section 9, the main points discussed in the review will be emphasized 
once again. 

2. The SU(N) 3D GG model 

The Euclidean action of the SU(2) 3D GG model reads [26] 



S= d^x 



1 '— ^2 1 .„ ^„.2 , A (t^a^2 2^2 



4,2 v^;.) + 2 (^.«&T + 4 (('I>T - ^^ 



(2.11 



where the Higgs field <I>" transforms by the adjoint representation, i.e. 
D^$" = d^<^'^ + e"'^'^A^^'^. The weak-coupling regime g^ <^ raw, which 
will be assumed henceforth, parallels the requirement that r/ should be large 
enough to ensure spontaneous symmetry breaking from SU(2) to U(l). At 
the perturbative level, the spectrum of the model in the Higgs phase consists 
of a massless photon, two heavy, charged VF-bosons with mass myY = 5^, 
and a neutral Higgs field with mass mn = f/\/2A. 

What is, however, more important is the nonperturbative content of the 
model, represented by the famous 't Hooft-Polyakov monopole [27]. It is a 
solution to the classical equations of motion, which has the following Higgs- 
and vector- field parts: 

• $" = (5"^ti(r), ii(0) = 0, u{r) — > t] — eiqy{—mHr) / {gr); 

Hf, = e^,yxd,yAl = -f - 4-it5{xi)6{x2)0{x3)S^3; 

• as well as the following action Sq = -^^. Here, k = g'^/mw is the 
weak-coupling parameter, e = e {niH/'mw) is a certain monotonic, 
slowly varying function, e > 1, e(0) = 1 (BPS-limit) [57], 

e(oo) ~ 1.787 [58]. 
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The approximate saddle-point solution (which becomes exact in the BPS- 
limit) was found in Ref. [26] to be 

q2 -^ / „ „ p-mH\Za-Zb\\ 

S = MSo + ^y T^^ - ^- ^ + 

8^ „^i \\Za-Zb\ \Za-Zb\ J 

where m^ -^ R -^ \za — zi)\, ggm = 47r, [gm] = [mass]^-*^'^. Therefore, while 
at niH -^ oo, the usual compact-QED action is recovered, in the BPS-limit 
one has 

8vr fj^^ \za - Zb\ 

i.e. the interaction of two monopoles doubles for opposite and vanishes for 
equal charges. 

When rriH < oo, the summation over the grand canonical ensemble of 
monopoles has been performed in Ref. [59] and reads 

Zmon = 1 + 5^ T7[ n / ^^^« Y. ""'^ = ^X^V'X 
\r -\ ' — 1 ^ .. LI J 



X exp \- \ d X 



J\f=l • a=l ■" qa=±l 

^ '(2.2) 

Here, x is the dual-photon field and ip is the field additional with respect to 
compact QED, which describes the Higgs boson. Furthermore, the monopole 
fugacity (i.e. the statistical weight of a single monopole), C, has the following 
form [26]: 

7/2 

Hi 
9 



^ = s!I^e-''>. (2.3) 



The function 5 = 6 {mu/myY) is determined by the loop corrections. It is 
known [60] that this function grows in the vicinity of the origin (i.e. in the 
EPS limit). However, the speed of this growth is such that it does not spoil 
the exponential smallness of C, in the weak-coupling regime under study. 

Our next aim will be to construct the SU(A^)-generalization of the par- 
tition function (2.2). Let us first introduce the {N — l)-dimensional vector 
H of the mutually commuting diagonal generators of the group SU(A^). 
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14 D. Antonov and M. C. Diamantini 

Together with certain pairwise hnear combinations of the off-diagonal gen- 
erators which, in analogy with the SU(2)-group, are called step (rising and 
lowering) generators E±i, i = 1, . . . , — 2~ ^ ^^^ diagonal generators form 
the following algebra: 



H,E±i 



(l±iE±i 



Ei , E^ 



m ■ 



Vectors Qi's here are called root vectors of the group SU(A^). The vector po- 
tential, ^^, a = 1, . . . , N'^ — 1, can be respectively decomposed into photons 
and W-bosons as 



^M = E[(^.^)^^- + (^/: 



+ A^H 



N(N~l)/2 

where from now on ^ = ^ . Next, embeddings of the SU(2)-monopole 

i i=l 

into the maximal Abelian subgroup [/(I) of the group SU(A^) can be 

characterized by the space-time variations of the Higgs field <J>". Outside 
the monopole core, $" can be chosen along the z-axis: <I>"(0, 0,z) — > fjaH. 
Here, fja is the (A^ — l)-dimensional vector of v.e.v.'s of the Higgs field $". 
For an arbitrary space-time direction, 

$"(r, 9, (^) = Q{e, (^)$"(0, 0, z)n-'^ {6, if). 

The matrix fi can be chosen such that 

<l>"(r,^,(p)=y« + Xf(r)7;-. (2.4) 

r 

Here, the 3-component matrix-valued vector 

rp; f Ei + E-i Ei — E-i ^ 

(where "i" under the square root denotes ^— 1) characterizes the embedding 
of the SU(2) Lie algebra into the root space of the group SU(iV), associated 
with the root gj. There exist therefore 2-^(-^— 1) embeddings, corresponding 
to the same number of monopoles. The constant Y"- = {fJa — {f]aQi)Qi) H, 
which parametrically depends on i, Y"" = y,"s, is the hypercharge associated 
with the i-th embedding. Since the vector r/a — {fjaqi)qi belongs to the plane 
containing qi and fja (and is perpendicular to gi), the first term on the r.h.s. 
of Eq. (2.4) breaks SU(iV) down to SU(2)xU(l)^-2, whereas the second 
term, with X^ — > Vaqi-, further breaks SU(2) down to U(l). Again, as in 
the SU(2)-case, the only part of the monopole vector potential which does 
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not vanish exponentially at large distances is the diagonal (photonic) one. 
It is given by 

g r 

that corresponds to the magnetic field 



HI 



9n 



Xi. 



.qM0,v)Hn-\e,ip), 



47rr^ 

where again g.^ = Air/g. Therefore, in the SU(A^)-case, monopoles also inter- 
act by means of the long-ranged Coulomb forces. This inter- monopole inter- 
action, mediated by dual photons, results in the following SU(iV)-analogue 
of the partition function (2.2): 



Z, 



N 



VxDil) exp 

2 

I 
2 



<rxx 



X ( \ {d^.xf + \ {d.i^f + ^V' - 2Ce<^-^ Yl ^°« (a^m 



(2.5) 



Here, the dual-photon field is described by the (A^— l)-dimensional vector x. 
Averaging in Eq. (2.5) over the Higgs field by means of the cumulant 
expansion one gets in the second order of this expansion [64] : 



"^mon 



/ ^Xexp - (fxi -{d^xf - 2^^ cos {gmqiX) 

L \ i . 

+2^^ / cf xcfy^ COS {gmqiXix))IC{x-y) COS (gmqjXiy)) 



1,0 



(2.6) 



In this equation, ^ = C 6xp 



2 ^niH [''^W , 



is the modified fugacity (which 



can be shown to remain exponentially small as long as the cumulant expan- 
sion is convergent), /C(x) = c9mDm.u{x) _ -^^ ^^^^^ jj^ ^--^ _ e~'"l^l/(47r|x|) is 
the 3D Yukawa propagator. The Debye mass of the dual photon, stemming 

from Eq. (2.6) by virtue of the formula ^ qfq^ = y<^°^, a,(3 = 1, . . . ,N —I, 

i 

reads 



rriD = grnV^C 



l + ^I 



N{N -1) 



where I = J d?xK.{x). At mu ~ 'niw-> I is given by [64] 



I 



5- exp 



^-■mH/mw 



(2.7) 



(21 
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16 D. Antonov and M. C. Diamantini 

The parameter of the cumulant expansion is O {^IN'^Y ^y virtue of 
Eqs. (2.3) and (2.8), one can readily see that the condition for this parameter 

to be (exponentially) small reads A^ < exp \-^ (e )]. Approximating e 

by its value at infinity, we find A^ < e^'^''^. Therefore, when one takes into 
account the propagation of the heavy Higgs boson, the necessary condition 
for the convergence of the cumulant expansion is that the number of colors 
may grow not arbitrarily fast, but should rather be bounded from above by 
some parameter, which is nevertheless exponentially large. ^ A similar anal- 
ysis can be performed in the BPS limit, mu <C g^ ■ There, one readily finds 
I ~ {gm/mnf, and 



^/A^2 oc iV^ exp 



47r . 
K V^~2 



Approximating e by its value at the origin, we see that the upper bound for 
A^ in this limit is smaller than in the vicinity of the compact-QED limit and 
reads N < e''/". 



3. String tension of the fiat Wilson loop in the fundamental 
representation 

At zero temperature, the 3D GG model is a clear analogy of the 2D XY- 
model in its continuum limit (see e.g. Refs. [26,62]). The analogy is due 
to the fact that vortices of the 2D XY-model correspond to monopoles of 
the 3D GG model, whereas spin waves correspond to free (non-dual) pho- 
tons. Disorder in both theories is produced by topological defects, i.e. by 
vortices or monopoles.-' Instead, spin waves or free photons cannot disorder 
correlation functions at large distances and affect them only at the distances 
smaller than the vacuum correlation length. As a result, monopoles lead to 
the area law of the Wilson loop, whereas free photons lead to its perimeter 
law. 

Thus, since confinement and string tension, we are interested in, are gen- 
erated by monopoles, photons will be omitted in this section. The partition 
function of the grand canonical ensemble of monopoles in the Euclidean 



' As for the convergence of tiie cumulant expansion itself, it is a natural requirement, which should 
be obeyed by any field theory with a normal stochastic, rather than the coherent, vacuum. 

■> At finite temperatures, disorder is primarily generated by W-bosons [52], which, at zero tempera- 
ture under discussion, are practically irrelevant due to their heaviness. The quantitative discussion 
of the role of W-bosons at finite temperature will be presented in subsection 8.4. 
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oo /-AT 

-^mon ~ 2-^ T7T \ ^'^^ 

X=0 ' 



-^ I d^xd^yff^{x)D,{x-y)(f^{y) 



(3.1) 



Here, Dq{x) = l/(47r|x|) is the 3D Coulomb propagator, and the monopole 

density is defined as p^{x) = Yl Qik^i^ ~ ^k) at AA > 1 and p^~^{x) = 0. 

fc=i 
Furthermore, the average is defined as 



(0)r 






n=0 



E 



O . 



(3.2) 



=±i,...,±j^a4^ 



Upon expUcit summation, the partition function (3.1) can be represented 
in the sine-Gordon-type form, which is the large-m/^ limit of the partition 
function (2.5): 



rN 



Pxexp 



d'^[\{d,xf 



2C ^ cos {gmqiX) 



(3.3) 



The Debye mass (2.7) becomes reduced to [64]: niD = gmV^C- 

The following comment is worth making at this point. The description 
of the grand canonical ensemble of monopoles in terms of the dual-photon 
field essentially implies the validity of the mean-field approximation. This 
approximation, which enables one to disregard fluctuations of fields of in- 
dividual monopoles, is only valid if the number of monopoles contained in 
the Debye volume, m^ , is large. Up to exponentially small corrections, 
the mean value of the monopole density, evaluated according to the for- 
mula pmean = yCT^ainT" ' ^^^^^ '^"^^^^ - C^(^ " 1), whcrC F^^) Jg ^^e 

three-volume occupied by the system. Therefore, the number of monopoles 
contained in the Debye volume is 



_3 



N 



1 



glVNC 



This is indeed an exponentially large quantity, even at iV '^ 1 . With the 
increase of N, the accuracy of the mean-field approximation is being further 
enhanced. 

Next, one can introduce the monopole field-strength tensor F^ which 
violates the Bianchi identity as \e^j_y\d^F^^ = gmf^ ■ Recalling that photons 
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are omitted throughout this section, we obtain 

Ffu (x) = -gms^uxdx I d^yDo {x - y) f {y) . (3.4) 

We will further use the formula 

N 



tr exp ( iOH j = N^ exp ( iOfla ) 



a=l 

where O is an arbitrary (A^ — l)-component vector and /2a's are the weight 
vectors of the fundamental representation of the group SU(iV), a = 1, . . . ,N. 
Together with the Stokes' theorem, this formula yields the following expres- 
sion for the Wilson loop defined at the AA-monopole configuration: 

^(C^)ln = ^ tr exp (!# f d^xFli^,^ = ^ £ W^.^ • (3.5) 

^ ^ a=l 

Here, 



W!f ^ exp (^|/2, j d^xFf^^^J^ 



(3.6) 



and S^jy (x) = / da^i, {x{S,)) 6 {x — x(^)) is the vorticity tensor current de- 

S(C) 
fined at a certain surface S(C) bounded by the contour C and parametrized 
by the vector x(^) with S, = (Co)Ci) standing for the 2D coordinate. A 
straightforward calculation with the use of Eq. (3.4) yields 

Wf = exp (ipa J d^xff^v) , (3.7) 

where rj (x; C) = f da a {x(S,)) du .^_\,^.. is the solid angle under which the 

surface S(C) is seen by an observer located at the point x, da^ = ^£^u\d(Ju\. 
(Note that, owing to the Gauss' law, rj = Air for a closed surface, i.e. when 
C is shrunk to a point.) The ratio of two W^^s defined at different surfaces, 
Si and S2, bounded by the contour C, reads 

^-^ ^^'^ - exp ( i^„ f da^^ixiO) [d'xp^{x)d, ^ '- 



^ S1UE2 ^ 

l[ew(-ijlaq^, I (ia^(x(0)5j«)^^^y (3.8) 



'^-'^ ^ S1US2 
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Due to the Gauss' law, the last integral in this equation is equal either to 
— 47r or to 0, depending on whether the point Zk lies inside or outside the 
volume bounded by the surface Si U S2. Since the product jlaQik is equal 
either to ±i or to 0, we conclude that pyjvL^| = 1. This fact proves the 

independence of W^ of the choice of the surface S in the definition (3.6). 
Clearly, this is the consequence of the quantization condition ggm = 47r used 
in the derivation of Eq. (3.7). 

A representation of the partition function (3.1), alternative to Eq. (3.3) 
and more appropriate for the investigation of the Wilson loop, is the one 
in terms of the dynamical monopole densities [64]. It can be obtained by 
multiplying Eq. (3.1) by the following unity: 



Vp5{p-p 



1^f\ 



VpDxew 



^g-n 



dPxx{p-f) 



so that the field x pla-ys the role of the Lagrange multiplier. We obtain for 
the partition function: 



Z, 



N 



VpDx^^V 



X 



N=0 
where the last sum is equal to 



^ / d^xcfypix)Do{x-y)p{y) + ig„, / d^xxp 



(3.9) 



exp 



(3.10) 



2C / (fx ^ cos igmqiX) 

i 

Accordingly, Eq. (3.7) becomes replaced by 

W^ ^Wa = -^ I VpDxexp l-^l d?xd'yp{x)D^{x - y)p{y) + 



and the full expression for the monopole contribution to the Wilson loop 

N 

[instead of Eq. (3.5)] reads l^(C)mon = ^ E ^a- 

a=l 

Let us next introduce the magnetic field according to the formulae 
d^B^ = p, e^^xduBx = 0. This yields 
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Wa = -^ J VB^5 (e^^xduBx) J Vxexp I J d^^[- ^Bf,+ 
+i9mXdtiBfj, + 2C ^ cos (gmqiX) + '^T^iUa I da^B^ > , (3.11) 

* E(C) ^ 

where, as has been discussed, the surface S(C) is arbitrary. The magnetic 
and the dual-photon fields can be integrated out of Eq. (3.11) by solving 
the respective saddle-point equations. From now on in this section, we will 
consider the contour C located in the (x, y)-plane. This naturally leads 
to the following Ansatze for the fields to be inserted into the saddle-point 
equations: B^ = 6f^3B{z), x = xi^)- For the points {x,y) lying inside the 
contour C, these equations then read: 

igmX + 9mB - Ampia5{z) = 0, (3.12) 

iB' -2C,^qi sin (gmqiX) = 0, (3.13) 

i 

where ' = d/dz. [Note that, differentiating Eq. (3.12) and substituting the 
result into Eq. (3.13), we obtain the equation 

x" - 25'mC XI ^* ^^"^ (.9mqiX) = gf^aS'iz), 

i 

which is the SU(A^)-generalization of the saddle-point equation obtained in 
Ref. [26].] These equations can be solved by using one more natural Ansatz 
for the saddle points of the fields, B{z) = flaB{z), xiz) = flaxi^)-, that, due 
to the formula 

/^a^fe = - Uab - -^ j , (3.14) 

makes Wa a-independent. Let us next multiply Eq. (3.13) by fla^ taking into 
account that, for any a, {N — 1) positive roots yield the scalar product with 
jla equal to 1/2, while the others are orthogonal to /!«. Equations (3.12), 
(3.13) then go over to 

2i(l)' + gl^B = A-Kidiz), B' + 2iCN siiKp = 0, (3.15) 

where (p = gmX/"^- One can readily check that the solution to this system of 
equations has the form 
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Biz) 



.SniD 



-mD\z\ 



Sm 



I _)_ Q-2m.o\z\ ' 



(t>{z) = 4 sgn z ■ arctan I e 



-mulzl 



(3.16) 



Taking the value -B(O) for the evaluation of the Wilson loop, we obtain for 
the string tension: 

N -1 AruD Stt N 



4tt ■ 



2N 



9n 



9ri 



'-V-C 



(3.17) 



This result demonstrates explicitly the nonanalytic dependence of a on g, 
but does not yield the correct overall numerical factor. That is because 
B{z), necessary for this calculation, is defined ambiguously. The ambiguity 
originates from the exponentially large thickness of the string: \z\ < d = 
m^ . It is therefore unclear which value of B{z) we should take: either 
B{0), as was done, or the one averaged over some range of z. It turns out 
that, in the weak-field (low-density) limit, the problem can be solved for 
an arbitrarily-shaped surface, by using the representation in terms of the 
Kalb-Ramond field, which will be described in the next section. 



4. SU(Ar) confining strings 

4.1. Fundamental representation 

The dual photon can alternatively be described in terms of the Kalb-Ramond 
field. One of the ways to introduce this field is to consider it as the field- 
strength tensor corresponding to the field S^, namely Bf^ = 2^£ij,i,\hv\- 
The Wilson loop (3.11) then takes the form'^ 



Wa 



1 



^£. 



Vh^y5 ( d^h^y 



Vxexp 



(Tx 



+ 2 X^t^uxdf^Kx + 2C X] ^°^ (.dmqiX) 



\Ku + 






WC^jy/J^jy 



S(C) 



(4.1) 

The field x can further be integrated out by solving the saddle-point equation 
of the form (3.13), where B' is replaced by 2^£ij,v\d^hy\ (and x depends 



^ Note the following correspondenee between tiic Coulomb interaction of monopole densities and 
the actions of the magnetic- and Kalb— Ramond fields: 

„2 r ^2 



4 7 ^"^2 7 ^2 



d-^xd^yp{x)Doix- y}p{y} 
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on all three coordinates). Using the Ansatz h^i, = ilahfj,u, we arrive at the 
following substitution in Eq. (4.1): 

/ d^x - xe^^xdiiKx + K^ cos (gmqix) =^ 

-2{N - 1)C / d^x fi^aarcsinhi^a - y/l + H^ . (4.2) 



Here, 



i^^u 



^" = 2ff„,C(iV-l) ^-^^-^^^^-^ ' (^-^^ 



and V [hf^u] is the multi-valued potential of the Kalb-Ramond field (or of 
monopole densities). Similarly to the case of compact QED [34], it is the 
summation over branches of this potential, which yields the summation over 
world sheets S(C) in Eq. (4.1). 

Another way to discuss the connection of the Kalb-Ramond field with the 
string world sheets is based on the semi-classical analysis of the saddle-point 
equations stemming from Eq. (4.1). These are analogous to Eqs. (3.15) and 
read: 

2id^cP + ^e^^xKx = 47riS^ , (4.4) 

7^—e^^uX^^Kx + 2iCN sin ,/. = . (4.5) 

Here, S^ = ^e^^xTji^x, and the field (j) is defined by the relation x = -^i^a'P- 
It is this auxiliary field 0, which establishes the correspondence between the 
Kalb-Ramond field and the stationary surface in the present approach. The 
summation over branches of the potential V\h^y\ is now replaced by the 
following procedure. One should restrict oneself to the domain |0| < vr and 
solve Eqs. (4.4), (4.5) with the conditions 

jim (l){x) = 0, (4.6) 

|x|— +00 

lim {(/. [m + en{i)\ - 4> \x{i) - en(C)]} = 27r, (4.7) 

e-+0 

where n(^) is the normal vector to S at an arbitrary point x(^). After that, 
one should get rid of the so-appearing S-dependence of Wa by extremizing 
the latter with respect to x(^). Owing to the formula 

(5 / day,yhy,y [f(^)] = / do ^yU ^yx \x{i)\ SxxiO , 
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such an extremization is equivalent to the definition of the extremal surface 
through the equation Hf^i,\ [x{^)] = 0. Due to Eqs. (4.5), (4.7), Hfj,,^\ indeed 
vanishes on both sides of the stationary surface. The reason for that is the 
following distinguished property of the extremal surface: (p is equal to vr on 
one of its sides and to — vr on the other, whereas for any other surface, \(j)\ is 
smaller than vr on one side and larger than vr on the other side. Therefore, 
according to Eq. (4.5), on both sides of the extremal surface, Ef^^xd^hux = 
0. This is equivalent to the equation e^jyA-f^/^i^A = and hence (upon the 
multiplication of both sides by e^'u'y) to H^^x = 0. [In particular, for a flat 
contour the extremal surface is the flat surface inside this contour. This can 
explicitly be seen from the solutions (3.16) of Eqs. (3.15). Namely, at any 
z such that \z\ <C d, 4){z) ~ vr • sgnz, B'{z) ~ 0, where both equalities hold 
with the exponential accuracy. The equation B'{z) ~ is equivalent to the 
condition of vanishing of H^i,x on both sides of the flat surface.] 

Next, the general form of an antisymmetric rank-2 tensor field h^^ is 

K^^ = d^A^ - d^Af, + e^^A^Af? . (4.8) 

The constraint d^hfj_y = 0, imposed in Eq. (4.1), is equivalent to setting ^^ 
to zero. From now on, we will promote /i^jy to include also the fields of 
free photons, ^^, by abolishing this constraint. Let us further go into the 
weak- field limit, \Ha\ <^ 1. Using the Cauchy inequality, we have 

' '•'-c(iv-i) cvmN - 1) ' ^^^ 

so that the weak-field limit is equivalent to the following low-density approx- 
imation: \p\ <^ (^^y2N{N — 1). To understand in which sense this inequality 
implies the low-density approximation, recall the mean monopole density we 
had from Eq. (3.3): 

|Plmean=^CA^(^-l). (4.10) 

Therefore, at large- iV, the low-density approximation implies that \p\ should 
be of the order N times smaller than its mean value. 

In the weak-field limit and with the constraint dfj_h^i, = removed, we 
now have the following expression for the total Wilson loop: 

X exp I - / d^x Y^^Hf^uX + l^u - ^l^ah^u^^u I • (4.11) 
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Here, ^*°* is given by the same integral over h^i,, but with S^,y set to zero, 
and -ff/ijyA = d^hy\ + d\h^^ + di,h\^ is the Kalb-Ramond field-strength ten- 
sor. [The apparent S-dependence of the r.h.s. of Eq. (4.11) is due to the fact 
that in course of the /i^j^-expansion of V \h^u\ , only one branch of this po- 
tential has been taken into account. As discussed above, the S-dependence 
disappears upon the summation over all the branches.] Owing to Eq. (4.^ 
Eq. (4.11) is obviously factorized as W {C,T.f^ll^_^^^^ = WaWl 
the free photonic contribution to the Wilson loop reads 



phot 



where 



^phot ^ g^p 



-^ '^ dXfj_ j) dx'^Do (x — x) 



AN 



c 



c 



(4.12) 



Doing the integration over h^i, in Eq. (4.11), we obtain 



w{c,^CL:-f^^M = ^-p{ -a' 



N - 1 



dx. 



C C 



dx D„ia{x - x) + 



+- 



m 



D 



dcJ^uim) / da^,{x{C'))Dmn {m - xiC')) 



S(C) 



S{C) 



(4.13) 



Note that the free photonic contribution completely cancels out of this ex- 
pression, i.e. it is only the dual photon (of the mass mu) which mediates 
the C X C- and S x S-interactions. 

One can further expand the nonlocal string effective action 



Oo- 



{gm 



D 



,N -I 

8N 



d(Tf,u{x{£,)) I daf,u{x{C))Drna {x{i) - x{i')) 



S(C) 



nc) 



(4.14) 

in powers of derivatives with respect to the world-sheet coordinates ^a^- 
Note that the actual parameter of this expansion is l/(m£)i?)^, where R ~ 
Y^Area(S) is the size of S (see the discussion below). The resulting quasi- 
local action is (cf. Refs. [36,65]): 

5str = cry d'eVS + «~' /^'eVgg"'(5aV)(56V) + 

+K j d^^^n + O f 3^75^1 . (4.15) 



^%. 



rWhR^ 



Here, da = d/d^"', and the following quantities characterize S: gabiO 
{daX^{C)){dbXf^{$,)) is the induced-metric tensor, g = det||g"^||, t^^{^) 
£"'^{daX^j,{Cj){dhXu{i))/ y/g is the extrinsic-curvature tensor, and TZ 
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(d"'da In y^) / yjg is the expression for the scalar curvature in the confor- 
mal gauge gab = ^/gSab■ 

The third term on the r.h.s. of Eq. (4.15) is known to be a full derivative, 
and therefore it does not actually contribute to the string effective action, 
while the second term describes the so-called rigidity (or stiffness) of the 
string [31,32]. The reason for the notation q~^, introduced in Ref. [31], is 
that, as has been shown in that paper, it is a which is asymptotically free. 
This asymptotic freedom then indicates that the rigidity term can only be 
infrared relevant if the respective /3-function has a zero in the infrared region. 
However, such a zero has not been found. Due to this fact, the rigidity term 
is not a good candidate to solve the old-standing problem of crumpling of 
large world sheets in Euclidean space-time. This necessitates seeking other 
possible solutions of this problem. One such solution, based on the string 
^-term, has been proposed in Ref. [31]. A possible derivation of such a term, 
within the SU(A^)-analogue of the 4D compact QED will be presented below. 

The string coupling constants in Eq. (4.15) read 

a = 271^^^-^, 4.16 

VN 9m 



A comment is in order regarding the negative sign of a. Up to a total 
derivative, the rigidity term reads 

where the Laplacian associated to the metric g"^ acts on a;^(^) as Ax^ = 
-j=da (\/gg SftX^). In the conformal gauge, one thus readily gets for the 

rigidity term J (P$^—r= (5^x^) . For a nearly flat surface, g ~ const = c^, 
and we have for the string propagator (in 2D Euclidean space): 

{xJo^um = ^-^ [d'p 



For Q < 0, this propagator has a tachyonic pole, whereas an unphysical mass 
pole shows up otherwise. We therefore conclude that the negative sign of a 
is important for the stability of strings. 

In the weak-field (or low-density) approximation, we can thus fix the pro- 
portionality factor 27r^ at a, that is close to the factor Svr of Eq. (3.17) which, 
however, could have not been fixed within the method of section 3. Another 
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important fact is that (in the weak-field approximation) the string tension 
and higher string couphng constants are the same for ah large enough sur- 
faces 5](C). [The words "large enough" here mean the validity of the inequal- 
ity R^ d, i.e. the string length R should significantly exceed the exponen- 
tially large string thickness d. Indeed, the general n-th term of the derivative 
(or curvature) expansion of S'str is of the order of am'j^R'^{R/d)^'^"^.] It is, 
however, worth noting that some of the terms in the expansion (4.15) van- 
ish at the surface of the minimal area corresponding to a given contour C 
(e.g. in the case of a flat surface considered in section 3). This concerns, for 
instance, the rigidity term.' 

Let us finally demonstrate how to introduce the Kalb-Ramond field in a 
way alternative to its definition via the -B^-field and incorporating automat- 
ically the free-photonic contribution to the Wilson loop. It is based on the 
direct combination of Eqs. (3.1) and (3.7), that yields 

In. \^"P [-^ _ 



Wa 



oo 

E 

Af=0 



d^xd^yp^ {x)Dq{x - y)p^ {y) + 



X exp< 



d^x 



0^,0 -d^T] 

9m 



inon 
2 



"^inon 



"^C^ COS {gmQi^) 



(4.18) 



i£a_, 



X + f^V- O^c can further use the following identity, which 



where 

shows explicitly how the Kalb-Ramond field unifies the monopole and the 

free-photonic contributions to the Wilson loop: 



exp 



-- / d^x[d^0-^d^r] 



Or, 



.N -I 

AN 



dx 



c 



c 



dx,,Dn(x — X ) 



Vhau exp 



.3^ / 1 D , ^ .^. Q r ^9 



d^ ( 7^^ + T^'f^fiuxd^Kx - —flah^u^^u ) I • (4.19) 



One can prove this by showing that both sides of this formula are equal to ' 



exp 



- / d^x{df,0 + gpa^^) 



'Indeed, one can prove the equality g''''(9ai/jiy){9i,tpi/) 



dadfjXfi — r^jj9c3::fi, and F^j^ is the Christoflfel symbol corresponding to the metric g"*". On the 



(V°'VaXf,){V''Vi,Xf,), where VaT>tx^, 
''. On 1 
together 



ab"<:-^IJ-' """ ^ ab 

other hand, the surface of the minimal area is defined by the equation T>°'T>aXi^{i;) 

with the respective boundary condition at the contour C. 

™ Obviously, in the noncompact case, when monopoles are disregarded and h^iy = df^A,^ — d^A 

the r.h.s. of Eq. (4.19) yields the free-photonic contribution to the Wilson loop, Eq. (4.12). 
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Inserting further Eq. (4.19) into Eq. (4.18), we obtain 



Wa ^ w{cy^ 



tot 



1 



Z^' 



,ot 



Dh^ijD^exp 



d^x{-hl,+ 



+ 2 0^i^uxd^Jiu\ - 2C XI ^°'^ {9mqi0) - —flji^v^ 



flU 



flU 



(4.20) 



Denoting x = ~0i we arrive back to Eq. (4.1) with Z^ 
constraint d^h^i, = removed, as it should be. 



N 
mon 



Z*°* and the 



4.2. Adjoint representation 

Let us now extend the ideas of the previous subsection to the case of the 
Wilson loop in the adjoint representation. The charges of quarks in this 
representation are distributed along the roots, so that in the formulae for the 
Wilson loop, Eqs. (3.6) and (3.7), fla should be replaced by qi. Noting that 
any root is a difference of two weights, we will henceforth in this subsection 
label roots by two indices running from 1 to N, e.g. Qab = fla — P-b -'^ 
Equation (3.14) then leads to the following formula: 



QabQcd — 2 ("^ac + ^bd — ^ad 



^bc) 



(4.21) 



according to which the product QabQcd may take the values 0, ±i, ±1. Owing 
to this fact, the ratio (3.8) is again equal to 1, i.e. for adjoint quarks, whose 
charge obeys the quantization condition gg^ = 47r, the Wilson loop is as 
surface independent, as it is for fundamental quarks. 

The adjoint-case version of Eq. (4.1) with the constraint 5^/i^jy = abol- 
ished [or of Eq. (4.20)] has the form: 



W{C)\ 



1 



2:tot 

N 



Vhf.^Vxexp 



d^x 



Iro 



K„+ 



c,d=l 



+ -^Qab 



^■•'^lly 



"'^fiu'T'fiu 



(4.22) 



S(C) 



The saddle-point equation, emerging from the above formula in the course 
of integration over x ; 

N 



-^^fivxd^ihyx 



9mC X] ^c<^ ^™ iamqcdx) , 



(4.23) 



c,d=l 



° In particular, this makes explicit the number of positive roots, 
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is to be solved in a way similar to the one of the previous subsection, namely 
by using the Ansatz h^i, = qabh^u, X = QabX- Multiplying both sides of 
Eq. (4.23) by Qab, we obtain: 



N 



^5mS 



^pivX 



diJiKx = ^ QahQcd sin {gmqabqcdx) ■ 



(4.24) 



cA=l 



Let us now compute the sum on the r.h.s. of this equation by virtue of 
Eq. (4.21). To this end, imagine the antisymmetric N x A^-matrix of roots 
Serf's. Obviously, for a fixed root qab-, there is one root qd (equal to qab) and 
a negative symmetric to it, whose scalar products with qab are equal to 1 
and —1. Our aim is to calculate the number n of roots, whose scalar product 
with qab is equal to zt^ .° According to Eq. (4.21), these roots belong either 
to the rows c = a, c = 5, or to the columns d = a, d = b, which in total 
contain AN — 4 elements. Among these, two are diagonal and other two 
yield the scalar product ±1. The remaining n = 4(A^ — 2) roots are precisely 
those which yield the scalar product ±2, so that this product equals ^ for 
2{N — 2) of them and — ^ for the other 2(A^ — 2). The sum on the r.h.s. of 
Eq. (4.24) thus takes the form 



sin(20) - sin(-20) + 2{N - 2) 



1 

— sine 

2 



1 . . 
-sm( 



2 [sin(20) + {N -2) sin < 



where the field (j) has been defined after Eq. (3.15). Let us now consider the 
limit A^ ^ 1, in which Eq. (4.24) takes the form sine/) = iHab, where [cf. 
Eq. (4.3)] 

Hab = TTT 7 Qab ^uuX^fiKx ■ 

The analogue of the potential (4.2) in the same limit then reads 



V 



flU 



4NC / d^x 



Hab arcsinh H, 



ab 



l + H± 



The string representation of the adjoint Wilson loop in the large- A^ limit is 
therefore given by Eq. (4.22) with the substitution 



d-^x 



N 



-,X£fj.uxd^^Kx + C X] ^°® {dmqcdX) 



c,d=l 



-V 



^^v 



° For all other N ~ N—2 — n roots, the scalar product vanishes, and so does the r.h.s. of Eq. (4.24). 
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(and the symbol "I?x" removed). 

We are finally interested in the adjoint-case counterparts of Eqs. (4.16), 
(4.17), one can obtain in the weak-field limit \Hab\ <C 1 .^ In this limit, 
formula (4.11) is reproduced, with the substitution /!„ -^ qab- As a conse- 
quence, the ratios of adjoint-case values of string couplings o"adj, oi^d\^ ^'^'^ 
«^adj to the respective fundamental-case values, Eqs. (4.16) and (4.17), are 
equal to jj^ — 2. In particular, for the string tensions this ratio coincides 
with the known leading large- iV QCD-result (see e.g. Ref. [66]). However, 
unlike the fundamental string, the adjoint one is unstable at large distances 
due to the production of P^+M^~-pairs. The string-breaking distance, i?c) 
can be estimated from the balance of the string free energy, (Jadj-Rc, and the 
mass of the produced pair, 2m\\/ [50]. By virtue of Eq. (4.16), one obtains 
Re oc d/n. Therefore, the breaking length of the adjoint string is in a factor 
K~^ larger than its thickness. However, since the very existence of the string 
implies that its length is much larger than its thickness [cf. the discussion in 
the paragraph preceding Eq. (4.18)], the question of (non-)existence of the 
adjoint string becomes purely numerical. 

4.3. k-strings 

The large- A^ ideas discussed in the Introduction have recently found a novel 
realization in the studies of the spectrum of /c-strings in SU(A^) gauge theo- 
ries. A /c-string is defined as the confining flux tube between sources in higher 
representations, carrying a charge k with respect to the center of the gauge 
group, Zjv, i.e. representations with nonvanishing A^-ality. These sources 
can be seen as the superposition of k fundamental charges, and charge con- 
jugation exchanges k- and [N — A:)-strings, so that non-trivial /c-strings exist 
only for A^ > 3 ; *^ their string tensions 0"^ can, and should, be used to con- 
strain mechanisms of confinement [69,70]. Results for the values of a^^ can 
be obtained by various approaches. Early results suggested the so-called 
"Casimir scaling" hypothesis for the ratio of string tensions [71]: 

R{KN)^^-^=^-^^^^C{KN), (4.25) 

(Ti iV — i 



P Note that, using the formula Hab = TjvF^atP *nd the Cauchy inequahty, one gets the following 

analogue of Eq. (4.9): \Hab\ < 2NC ' This clearly leads to the same definition of the weak- field limit 
in terms of the low-density approximation, as in the fundamental case. Namely, the weak-field 
limit corresponds to densities \p\, which are of the order N times smaller than the mean one (4.10). 

'^ We will not consider here high-dimensional representations that are screened by gluons and do 
not yield a genuine asymptotic string tension. 
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where cJi is the fundamental string tension. These are based on the argu- 
ments of the Parisi-Sourlas dimensional reduction of the 4D QCD to the 2D 
one, which is due to stochastic vacuum fields, and the further use of the 
fact that, in 2D, confinement is produced by the one-gluon exchange. Re- 
cent studies in supersymmetric Yang-Mills theories and M-theory suggest 
instead a "Sine scaling" formula: 

R{k,N)= . \ ' . 4.26 

^ ^ sm(7r/iV) ^ ^ 

Corrections are expected to both formulae, but the form of such corrections 
is unknown for the physically relevant case of a 4D, nonsupersymmetric, 
SU(A^) gauge theory. 

In the large- A^ limit, where the interactions between flux tubes are sup- 
pressed by powers of 1/A^, the lowest-energy state of the system should be 
made of k fundamental flux tubes connecting the sources, hence: 

fc — fixed 

R{k,N) ""-^^ k . (4.27) 

Both the Casimir and the Sine scaling formulae satisfy this constraint; they 
also remain invariant under the replacement k ^^ {N — k), which corresponds 
to the exchange of quarks with antiquarks. However, it has been argued in 
Refs. [72, 73] that the correction to the large-A^ behavior should occur as 
a power series in 1/N'^ rather than 1/A^ . "^ Clearly, such a behavior would 
exclude Casimir scaling as an exact description of the /c-string spectrum. 

Recent lattice calculations have provided new results for the spectrum 
of /c-strings both in three and four dimensions [74-79]. They all confirm 
that Casimir scaling is a good approximation to the Yang-Mills results. To 
be more quantitative, one could say that all lattice results are within 10% 
of the Casimir-scaling prediction, and that deviations from it are larger in 
four than they are in three dimensions, in agreement with strong-coupling 
predictions [78] . The taming of systematic errors is a crucial matter for such 
lattice calculations, and it can only be achieved by an intensive numerical 
analysis. In four dimensions, the higher statistics simulations presented in 
Ref. [78] show that corrections to the Casimir-scaling formula are statistically 
significant, and actually favor the Sine scaling. Finally, it has been pointed 
out in Ref. [79] that higher-dimensional representations with common A^- 
ality do yield the same string tension, as expected because of gluon screening. 



" One should study with some care whether the arguments presented in Refs. [72, 73] hold inde- 
pendently of the space-time dimensionality. 
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These numerical results trigger a few comments on Casimir scaling. The 
original argument [71] was based on the idea that a 4D gauge theory in a ran- 
dom magnetic field could be described by a 2D theory without such a field. 
Besides the numerical results, there is little support for such an argument 
in QCD; moreover it is not clear that the same hypothesis could explain the 
approximate Casimir scaling observed in three dimensions. On the other 
hand, Casimir scaling appears "naturally" as the lowest-order result, both 
at strong-coupling in the case of fc-strings in the Hamiltonian formulation 
of gauge theories, and in the case of the spectrum of bound states in chiral 
models. Corrections can be computed in the strong-coupling formulation and 
they turn out to be oc {D — 2)/N - see e.g. Ref. [78] for a summary of results 
and references. While strong-coupling calculations are not directly relevant 
to describe the physics of the continuum theory, it is nonetheless instructive 
to have some quantitative analytic control within that framework. Last but 
not least, Casimir scaling also appears at the lowest order in the stochastic 
model of the QCD vacuum [80]. In view of these considerations, it is fair 
to say that approximate Casimir scaling should be a prerequisite for any 
model of confinement, that corrections should be expected, and that these 
corrections are liable to yield further information about the nonperturbative 
dynamics of strong interactions. Moreover, it would be very interesting to 
improve our understanding of some other aspects of the fc-string spectrum, 
like e.g. the origin of the Sine scaling for nonsupersymmetric theories, or 
the structure of the corrections to this scaling form. 

We will prove below that Casimir scaling occurs to a high accuracy in the 
low-density approximation of the SU(iV) 3D GG model [81]. The /c-string 
tension is defined by means of the k-th. power of the fundamental Wilson 
loop. The surface-dependent part of the latter can be written, in terms of 
the dual-photon field, as follows: 

TV / 

(^^(C))mon= E («^P 
ai,...,afe=l\ 

In this equation, 

S(C) 

where again S(C) is an arbitrary surface bounded by the contour C and 
parametrized by the vector x(^). 

The independence of Eq. (4.28) of the choice of S(C) can readily be seen 
in the same way as for the {k = l)-case. The S-dependence rather appears in 




d xT, {x) X {^) 



(4.28) 
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the weak-field, or low-density, approximation, which is equivalent to keeping 
only the quadratic term in the expansion of the cosine in Eq. (3.3). As has 
been demonstrated in subsection 4.1 [cf. the discussion around Eqs. (4.9), 
(4.10)], the notion "low-density" implies that the typical monopole density 
is related to the mean one, pmcmi = C^{^ ~ ^)i by the following sequence 
of inequalities: 

Aypieal « C " O(A^) « Pinean = ('0 (N^) . (4.29) 

We will discuss in some more detail the correspondence between the low- 
density and the large- A^ approximations after Eq. (4.36). 

Denoting for brevity a = —g J d^xT, [x)x{x), we can rewrite Eq. (4.28) 
as j^ 

{WkiC))^^^= J2 (e-('^-+-+^^"J), (4.30) 

( ai,...,a;,-l \ 

Vwith possible coincidences/ 

where in the low-density approximation the average is defined with respect 
to the action 



d'^x 



2 



1 /o ^\2 , "I'D ^2 



(4.31) 



Similarly to the fundamental representation, in this approximation the string 
tension for a given k is the same for all surfaces S(C), which are large enough 
in the sense \/S S> d, where S is the area of S(C). In particular, the fun- 
damental string tension found above reads o"i = -t^cj, where a = Air"^ ^^ , 
and the factor ^^ is the square of a weight vector. 

To evaluate Eq. (4.30) for A; > 1, we should calculate the expressions of 
the form 



(4.32) 



where [k — n) weight vectors fla 's are mutually different and also different 
from the vector /Iq. . By virtue of the formula flafib = 2 {^a^ ~ Iv) ' ^^ obtain 
for Eq. (4.32): 




^-1 / 2 , N 1 
[n +k — n) 

2N ^ ' 2N 



k—n—l 



2n{k-n) + 2 V I 



1=1 



k(N-k) 1 , ^ . 



2N 2 

(4.33) 

We should further calculate the number of times a term with a given n 

appears in the sum (4.32). In what follows, we will consider the case k < N, 

~m -- k\ 
''k n!(fc— n)! 



although k may be of the order of N . Then, Cu = f,t^ w possibilities exist 
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to choose out of k weight vectors n coinciding ones, whose index may acquire 
any values from 1 to A^. The index of any weight vector out of the remaining 
{k — n) ones may then acquire only (iV — 1) values, and so on. Finally, the 
index of the last weight vector may acquire (A'^ — k + n) values. Therefore, 
the desired number of times a term with a given n appears in the sum (4.32) 
is 

CJ^N ■ {k - n){N - 1) ■ {k - n - 1){N - 2) ■ ■ ■ 1{N - k + n) = 



k\N\ 
n\{n + N -k-l)l ' 



Cj:A%--+\k -n)\ = ^„„ , ,; •, ,„ , (4.34) 



where yl^~"+^ = ,jy_^^'^_^s, . Equations (4.33) and (4.34) together yield for 
the monopole contribution to the Wilson loop, Eq. (4.30): 






1 "2^ 



m(C)}„, = km<^-'^ . 2 _,_ e--^'^ . (4.35) 

n=l 

where C = -^m i^ proportional to the Casimir of the rank-A; antisymmetric 
representation of SU(iV). We have thus arrived at a Feynman-Kac~type 
formula, where, in the asymptotic regime of interest, S -^ cxd, only the 
first term in the sum is essential. The /c-string tension is therefore ak = 
Ca, yielding the Casimir-scaling law (4.25). It is interesting to note that 
the Casimir of the original unbroken SU(iV) group is recovered. This is 
a consequence of the Dirac quantization condition [61], which distributes 
the quark charges along the weights of the fundamental representation and 
the monopole ones along the roots. The orthonormality of the roots then 
yields the action (4.31), which is diagonal in the dual magnetic variables; 
the sum of the weights squared is responsible for the Casimir factor, since 
C = {flai + • • • + /Xflj. ) , where all k weight vectors are different from each 
other. Therefore, terms where all k weight vectors are mutually different 
yield the dominant contribution to the sum (4.30). Their number in the 
sum is equal to T^yirpry , that corresponds to the (n = l)-term in Eq. (4.35). 
Let us further address the leading correction to Casimir scaling, which 
originates from the non-diluteness of plasma. Expanding the cosine up to 
the quartic term and using the formula 

NiN-l)/2 ^ 



E 



Qr'i^<lIQ! = 2(N + l) i^""^'' + '^"'''^''' + '^"''^''^) 



^ 2(iV + l 
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which stems from the orthonormahty of roots, we obtain the action 



(4.36) 



By virtue of this formula, one can analyze the correspondence between the 
1/A^-expansion and corrections to the low-density approximation. The nat- 
ural choice for defining the behavior of the electric coupling constant in 
the large- A^ limit is the QCD-inspired one, g = 0{N^^''^). To make some 
estimates, let us use an obvious argument that the X'^^ld configuration 
dominating in the partition function is the one where every term in the ac- 
tion (4.36) is of the order of unity. When applied to the kinetic term, this 
demand tells us that the characteristic wavelength / of the field x is related 
to the amplitude of this field as / ~ |xl^^- Substituting further this estimate 
into the condition l^Tn'j-,\x\'^ ~ 1, we get |xP ~ "t-d- The ratio of the quartic 
and mass terms, being of the order of IxPffm/-^) can then be estimated as 
"'^"' = gfriV N ~ ^V(- With the exponentially high accuracy, this ratio 
is small, provided N <0 (e'^"'^).^ Therefore, the non-diluteness corrections 
are suppressed only for A^'s bounded from above by a certain parameter. 
However, due to the exponential largeness of this parameter, the obtained 
constraint leaves enough space for N to be sufficiently large to provide the 
validity of the inequalities (4.29). 

To proceed with the study of the non-diluteness correction, one needs to 
solve iteratively the saddle-point equation, which corresponds to the aver- 
age (4.30) taken with respect to the approximate action (4.36). Since it has 
been demonstrated above that the string tension is defined by the averages 
/g««(/^oiH \-^J.a^)\ -^j^ere all k weight vectors are mutually different, let us 
restrict ourselves to such terms in the sum (4.30) only. Solving then the 
saddle-point equation with the Ansatz X = Xo + Xi; where |xi| ^ |xo|) we 
obtain for such a term: 

2 I- 

- In /e*<'^-i+-+'^«fc)\ = ^C SxSyll {x) D^^ (x - y) S {y) + AcS . 

(4.37) 
The first term on the r.h.s. of Eq. (4.37) yields the string tension a^ = Ca, 



" This eonstraint is similar to those which were derived at the end of sect. 2 as the necessary 
conditions for the stochasticity of the Higgs vacuum. We therefore conclude that the conditions of 
stochasticity of the Higgs vacuum and of the validity of the dilute-plasma approximation parallel 
each other. 
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while the second term yields the desired correction. This term reads 

= y ^^^^/^^M- (^i) d^f^- (^2) daxp (X3) dax, (x,) d^' 8%^ df df^^H , 

(4.38) 

4 

where I = f d^x Y\ ^mu {^ ~ ^i)- The action (4.38) can be represented in 

1=1 
the form 



A5 = 5 / da^,y{xi)da^,y{x2)D{xi - X2) x 



X / da\p{x3)daxp{x4)D{x3 - X4) x G{xi - X3) . (4.39) 

Here, D and G are some positive functions, which depend on molxi — Xj\ 
and vanish exponentially at the distances > d. They can be represented as 
D{x) = mjyT){mD\x\), G{x) = m'\^Q{m]:)\x\), where the functions T) and Q 
are dimensionless.* 

The derivative expansion yields the Nambu-Goto action as the leading 
term: 



da puixi) da pu{x2)D{xi - X2) 



CTD j d\.M^) + o(^Y (4.40) 



Here, a^ = 2m?jj J d'^ z'D{\z\) (with z being dimensionless) . Recalling that 
daapix) = \/g{x)taf3{x)d'^^, we may further take into account that we are 
interested in the leading term of the derivative expansion of the action AS, 
which corresponds to the so short distance |xi — X3I, that t^^ (xi) tap (^3) ~ 
2. (Higher terms of the derivative expansion contain derivatives of t^^ and 
do not contribute to the string tension.) This yields for the integral in 
Eq. (4.39): . 

al J d\d\'.Jg{x^)g{x.,)G{x^ - X3) ~ 

~ ^jdaap (x 1) daap {x^) G (xi - ^3) = ^ ^g j d^iVi + O (^ 
where ac = 2m'j^ J d^zQ^zl). We finally obtain from Eq. (4.39) that 



* Our investigations can readily be translated to the stochastic vacuum model of QCD [80] for the 
evaluation of a correction to the string tension, produced by the four-point irreducible average 
of field strengths. In that case, the functions D and G would be proportional to the gluonic 
condensate. 
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^ {gCfalcjG _ a {gCfmp _ aC 

^^k — o/AT , 1N '^ A AT- , 1 ^~ 



2{N + l)m% 4 A^ + 1 iV + 1' 

where a is some dimensionless positive constant. This yields 

aC 



Gk + AcTfe = <tC 1 + 



iV + 1 



In the hmit A; ~ A^ ^ 1 of interest, the obtained correction is 0(1), while 
for k = 0{\) it is 0{1/N). The latter fact enables one to write down the 
following final result for the leading correction to Eq. (4.25) due to the non- 
diluteness of monopole plasma: 



R{k, N) + AR{k, N) = '^''^ '^'^ = c{k, N) 

CTi + A(Ti 



(fc-l)(iV-fc-l) 

2N{N + 1) 



This expression is invariant under the replacement k <-^ {N — k) just as the 
expression (4.25), which does not account for non-diluteness. The fact that, 
at /c ~ A^ ^ 1, the obtained correction to the Casimir-scaling law is 0(1), 
indicates that non-diluteness effects can significantly distort the Casimir- 
scaling behavior. 

5. Generalization to the SU(Ar)-analogue of 4D compact 
QED with the 0-terni 

In this section, we will consider the 4D-case and introduce the field- 
theoretical 0-term.^ As was first found for compact QED in Refs. [33, 36] 
by means of the derivative expansion of the resulting nonlocal string effec- 
tive action, this term generates the string ^-term. Being proportional to the 
number of self-intersections of the world sheet, the latter might be important 
for the solution of the problem of crumpling of large world sheets [3, 31]. In 
this section, we will perform the respective analysis for the general SU(A^)- 
case under study, in the fundamental and in the adjoint representations, 
as well as for /c-strings. It is worth noting that, in the lattice 4D compact 
QED, confinement holds only at strong coupling. On the other hand, the 
continuum counterpart [33, 36] [of the SU(A^)-version] of this model, which 
we are going to explore, possesses confinement at arbitrary values of cou- 
pling. However, we will see that the solution to the problem of crumpling 



" In the case when the ensemble of Abehan-projected monopoles is modeled by the magnetically- 
charged dual Higgs field, the effects produced by this term on the respective string effective action 
have been studied in Ref. [67]. 
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due to the O-term. is only possible in the strong-coupling regime, implied in 
a certain sense. Note also that the continuum sine-Gordon theory of the 
dual-photon field possesses an ultraviolet cutoff, A, which appears in course 
of the path-integral average over the shapes of monopole loops. However, 
the ^-parameter is dimensionless, and consequently its values, at which one 
may expect the disappearance of crumpling, will be cutoff-independent. 

The full partition function, including the ^-term and the average over 
free photons, reads [cf. Eq. (3.1)] 






where 
S 






i fd^xd'yj^ix)Do{x - y)j^iy) + ^ / d^x\j^, (5.2) 



and F^i, = dfj_Ay — d^A^. For AA = 0, the monopole current 
j^ is equal to zero, whereas for AA > 1, it is defined as j^ = 

N 
9m X^ Qik § dz^{T)5 (x — x^(r)). The couplings g and Qm are now dimen- 

fc=i 
sionless and obey the same condition gg^ = 47r as in the 3D-case. We 

have also parametrized the trajectory of the fc-th monopole by the vector 

1 
^^^l{'T) = y^ + -2^^(''")) where y^ = j dTX^^{T) is the position of the trajectory, 



whereas the vector zf^{T) describes its shape, both of which should be aver- 
aged over. The fugacity of a single-monopole loop, (", entering Eq. (5.1), has 
the dimensionality [mass]*^, C, oc e~ ™°", where the action of a single /c-th loop, 

obeying the estimate S'mon ot -^ f dry {z'^) , is assumed to be of the same 

'^ 
order of magnitude for all loops. Finally, in Eq. (5.1), Dq{x) = l/(47r^x^) is 

the 4D Coulomb propagator, and the average over monopole loops is defined 

similarly to Eq. (3.2) as follows: 



^■. — n "^ ^^(• AT 1 \ 



«=o-' i„=±i, ^m^^ 



The particular form of the path-integral average over the shapes of the loops, 
2;„(r)'s, here is immaterial for the final (ultraviolet-cutoff dependent) expres- 
sion for the partition function [see e.g. Ref. [68] for a similar situation in the 
plasma of closed dual strings in the Abelian-Higgs-type models]. The only 
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thing which matters is the normahzation (1)^ ,-. = 1, that will be implied 
henceforth. The analogue of the partition function (3.9), (3.10) then reads 



Z, 



N 



A„ 



/ ^^3^,ex];)^-S j^,A^ | 



X / ^X^exp 



d^x 



2CE 



cos I ^ ) + tXt^Jfi 



(5.3) 



with the action S given by Eq. (5.2), while the full partition function is 



Z 



N 



VA e--J '^^^^i'^ Z^ 

i^^^e "i f- ^mon 



A„ 



In Eq. (5.3), 



and A ~ \y^\ / [z^) > jz^p 



xlxl, ■ ■ ■ , \/ xf! ^X^ ^ 



(5.4) 



is the ultraviolet cutoff. Clearly, unlike the 
3D-case without the ^-term, the ^^-field is now coupled to j^, making Z^^^ 
^^-dependent. This eventually leads to the change of the mass of the dual 
photon Xfj.! as well as to the appearance of the string 9-term. 

Let us address the fundamental case first. After the saddle-point integra- 
tion over Xfj_ ,"^ the analogue of Eqs. (4.1), (4.2) takes the form 



Wa = ^ I T^h^u exp 



S 



fJ,U 






dcrfj.uhf_,u \ , (5.5) 



S(C) 



where the constraint dfj_h^i, 
Ramond action reads 



was already abolished. Here, the Kalb- 



''flU 



"^ 4 ^'^ 327r2 



h^,uh^,u + V 



''fJLV 



(5.6) 



The potential V in this equation is given by Eq. (4.2) with the symbol J d^i 
replaced by J d^x, and 



^ Unlike the superrenormalizable 3D-case, where loop correetions to the saddle point are rapidly 
converging, this is not necessarily the case in four dimensions. This fact is, however, clearly 
unimportant in the weak-field limit [Eq. (5.7) below], where the cosine in Eq. (5.3) is approximated 
by the quadratic term only, and the saddle-point integration over Xm becomes Gaussian. As we 
have seen in the previous section, this limit is already enough for the discussion of the string 
representation. One of the authors (D.A.) is grateful to H. Gies and E. Vicari for drawing his 
attention to this issue. 
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Ha 



5A 



C(iv-i) 



fJ-a 



Cf^i^^v 



In these formulae, O^u = 2^fj.iyXp^Xp} a-iid the absolute value is defined in 
the same way as in Eq. (5.4), i.e. again with respect to the Lorentz indices 
only. Note that the form to which the ^-term has been transformed is quite 
natural, since the respective initial expression of Eq. (5.2) can be rewritten 
modulo full derivatives as 



iOg^ 



(fxA^j 



•V 



87r2 
where [cf. Eq. (3.4)] 



327r2 



d''x{F,, + F^, 



F +F^ 



FffAx) 



e^uXpdx / d yDo{x - y)fp (y), d^F^^^ = ji . 



Next, the mass of the Kalb-Ramond field, equal to the Debye mass 
of the dual photon, which follows from the action (5.6), reads m£) = 



9V, 

An 



+ (2^) , where ij = yJNQ/ K. In the extreme strong-coupling 

limit, g -^ cxd, this expression demonstrates the important difference of the 
case = from the case 7^ 0. Namely, since r]{g) oc e~'^°°^*'^ -^ 1, 
rriD ^ at = 0, whereas m£) ^ cx) at 7^ 0. In another words, in the 
extreme strong-coupling limit, the correlation length of the vacuum, equal 
to d, goes large (small) at = (0 / 0). 

In the weak-field limit, \Ha\ <. 1, Eq. (5.5) yields [cf. Eq. (4.11)]: 

1 



X exp 



d^x 



2 



tot, a 
weak— field 



g 



12r/' 



t?2 

2 iJ.y\ 



Ap 



2:tot 
^2 



Vhf.uX 



iOg\ , 



ig 



l^ah^u^ 



p,u^p,u 



(5.7) 



where again Z is the same integral over h^u-, but with S^j^ set to zero. 
Integrating over h^y, we then obtain: 



W{C,^f:li:__^.,^ = e^V 



N -I 



dx. 



"'■^ H^rnjj \^ 



X') + 



c 



c 



+— I d xd x'Drno{x - x') S^j.(x)S^,,(x') + -— ^i;^,,(x)i;^j,(a;') 



Svr^ 



(5.8) 
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where D„i{x) = mKi{m\x\) /{Ait'^\x\) is the 4D Yukawa propagator with 
Ki denoting the modified Bessel function. [Note that, at ^ = 0, Eq. (5.8) 
takes the form of Eq. (4.13).] Further curvature expansion of the S 



IJLU 



T,f^i, interaction [analogous to the expansion of the S^j^ x S^jy interaction 
in the action (4.14)] yields the string 0-term equal to icfund'^- Here, u = 
^ J d'^^\'^E,"'^{9atfj.i/){dbi^u) is the number of self-intersections of the world 
sheet, and the coupling constant Cfund reads 

{N - 1)6 



Cfund 



8iV 



(5.9) 



As has already been discussed, Cfund is A-independent, since (similarly to 
the rigidity coupling constant a) it is dimensionless. We therefore see that, 
at 

N 



a fund 



1 



2N 



'N-V 

2N 



(5.10) 



Cfund becomes equal to vr, and self-intersections are weighted in the string 
partition function with the factor (— I)'', that might cure the problem of 
crumpling. This is only possible at g > g^^^'^ = 4. /-^JL_^w paralleling the 
strong-coupling regime, which should hold in the lattice version of the model 
(cf. the discussion in the beginning of this section). In the extreme strong- 
coupling limit, understood in the sense g ^ ^^fund^ only one critical value, 

whereas 0^^'^'^ -^ 0, i.e. 05^°'^ becomes a 
0. 



3 fund 



vr 



N-l 
2N ' 



e=o 



^fund^ survives, 

spurious solution, since Cfund 

With the use of the results of subsection 4.2, the adjoint-case large-A^ 
counterpart of Eq. (5.10) can readily be found. Indeed, Eq. (5.9) in that 
case becomes replaced by 

Cadj 



+ (^J 



so that Cadj is equal to vr at 



^adj 



vr 
2 



1± 



m 



adj 



Note that the respective lower bound for the critical value of g, (77v>i ~ 
4^/^: ,^ is only slightly different from the value 57v>i ~ 4v2vr. Similarly 



nfund 



" Clearly, at g = g'" 

^ As in tile previous footnote, at tlie particular value of <;, 3 = g^L^^, d'^^ = 9^ 



N-l 

4N 



adj 



/jadj 



nadj 
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to the fundamental case, at g ^ 9n^i^ ^- ^ becomes a spurious solution, 
whereas 6^^ ^ tt. Thus, at iV S> 1 and in the strong-coupling limit, un- 
derstood in the sense g ^ 5'7v>i' ^^^ critical fundamental- and adjoint-case 
values of 9, at which the problem of crumpling might be solved, are given 
by the following simple formula: 0^^^'^ = 2^+'' = §• 

Finally, for fe-strings, the ratio Cfc/cfund is the same Casimir one, C{k, N), 
as the ratio of string tensions, as long as A^ < O (^e"^™""'^).^ This enables 
one to readily find 9^ as well. Namely, 



5fc _ ^ 



N - 1 



f N-l Y /167r 



2NC{k,N) V \2NC{k,N)) \ g^ 



that is only possible at g > gk ^ g^"^""^ ^J C {k , N) > g^"^""^. 

6. Geometric aspects of confining strings: the physics of 
negative stiffness 

The role of antisymmetric tensor field theories for confinement has been 
studied in detail in Ref. [33]. By analyzing compact antisymmetric tensor 
field theories of rank h — 1 Quevedo and Trugenberger have shown that, 
starting from a "Coulomb" phase, the condensation of (d — h — l)-branes 
(where D = d + 1 is the space-time dimension) leads to a generalized con- 
finement phase for {h — l)-branes. Each phase in the model has two dual 
descriptions in terms of antisymmetric tensor of different ranks, massless 
for the Coulomb phase and massive for the confinement phase. Upon inte- 
gration over the massive antisymmetric tensor field (in the case of an even 
number of dimensions - when the ^-term is absent), one obtains the string 
effective action (4.14). The derivative expansion of this nonlocal string effec- 
tive action, up to the term next to the rigidity one, produces the following 
confining-string action: 



S = I d'e^/g g^'VaX^ (t - sV^ + -^vA V,x^. 



(6.1) 



Here, Va is the covariant derivative with respect to the induced metric gab 
idaXfj.){dhX^) on the surface parametrized by x(^O)Ci)- 



^ In 4D, the ratio of string tensions indeed equals C{k, N) as in 3D, as long as the dual-photon 
field can be treated as a free massive field, since the rest of the derivation of Eq. (4.35) is based 
on the properties of fundamental weights only. 
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In (6.1), the first term provides a bare surface tension 2T, while the second 
accounts for rigidity with stiffness parameter s. The last term can be written 
(up to surface terms) as a combination of the fourth power and the square of 
the gradient of the extrinsic curvature matrices, with M being a new mass 
scale. It thus suppresses world-sheet configurations with rapidly changing 
extrinsic curvature; due to its presence, the stiffness s may not necessarily 
be negative, as in Eq. (4.17) above, but also positive. We analyze the model 
(6.1) in the large-I? approximation. To this end, we introduce a Lagrange 
multiplier matrix L"* to impose the constraint gab = {daX^){di)X^), extending 
the action (6.1) to 

S + j Cpi^ L^'idaX^dbX^-gab) . (6.2) 

Then we parametrize the world sheet in a Gauss' map by x^{S,) = 
(^o,ei,</''(0), ii = 2,...,D-l), where -/3/2 < ^o < /3/2, -R/2 < ^ < 
R/2, and (j)^{^) describe the D — 2 transverse fluctuations. With the usual 
homogeneity and isotropy Ansatz gab = pSab , L"-^ = Lg"-^ of infinite surfaces 
(/9, R -^ oo) at the saddle point, we obtain 

S =2 I Si[T + L(l -p)] + I cPi daCp'V (T, s, M, L, V^) dacP' , (6.3) 

where 

V{T,s,M,L,V^)=T + L^sV'^ + ^V\ (6.4) 

Integrating over the transverse fluctuations, in the infinite-area limit, we get 
the effective action 

5^^ = 2Aext[r + L(l -/>)] + 

+Ae.t^^^P [(fpin[p''ViT,s,M,L,p)] , (6.5) 

where ylext = PR is the extrinsic, physical, space-time area. For large D, 
the fluctuations of L and p are suppressed and these variables take their 
"classical values" , determined by the two saddle-point equations 

= f{T,s,M,L) , p= \. n ' (6-^) 

f{T,s,M,L) 

where the prime denotes a derivative with respect to L and the "saddle- 
function" / is defined by 

/ (T, s,M,L) = L- ^^ J dp pin [p^V (T, s, M, L, p)] . (6.7) 
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Using (6.6) in (6.5) we get 5^^ = 2 (T + L) ^ext showing that T = 2{T + L) 
is the physical string tension. 

The stabihty condition for the Euchdean surfaces is that V (T, s, M, L,p) 
be positive for all p'^ > 0. However, we will require the same condition 
also for p"^ < 0, so that strings propagating in Minkowski space-time are 
not affected by the propagating states of negative norm which plague rigid 
strings. The stability condition becomes thus \/T + L > \sM/2\, which 
allows us to introduce the real variables R and / defined by 

In terms of these, the kernel V can be written as 

M^V{T,s,M,L,p) = {R^ + ff + 2{R^-I^)p^+p'^ . (6.9) 

In order to analyze the geometric properties of the string model (6.1) 
we will study two correlation functions. First, we consider the orientational 
correlation function 

9ab{i-O^{da^\0db<t>\e)) 

for the normal components of tangent vectors to the world sheet. 

Secondly, we compute the scaling law of the distance dE in embedding 
space between two points on the world sheet when changing its projection d 
on the reference plane. The exact relation between the two lengths is 

dl = d^ + Y.{W{i)-(t>M')\'')- (6.10) 

i 

When d\ oc d^ this implies that the Hausdorff dimension of the surfaces, 
defined as d\ = (d^)^'^^, is equal to 2, and the surface is smooth. 

In order to establish the properties of our model, we analyze the saddle- 
point function f {T,s,M,L) in (6.7). To this end, we must prescribe a 
regularization for the ultraviolet divergent integral. We use dimensional 
regularization, computing the integral in (2 — e) dimensions. For small e, 
this leads to 

/(r,,,M,L) = L + ^(fl^-/^)ln:?!±^-^fl/(| + aaCa.4^ 

(6.11) 

where A = /x exp(2/e) and ;U is a reference scale which must be introduced 
for dimensional reasons. The scale A plays the role of an ultraviolet cutoff, 
diverging for e ^ 0. 
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The saddle-point function above is best studied by introducing the di- 
mensionless coupHngs t = T/A^, m = M/A and I = L/A^. We will study 
in detail the case s = 0, in which the saddle-point equations can be solved 
analytically since R = I. This choice is not too restrictive since, as we will 
show, s = is the infrared fixed point. We get 



where c = l/(327r^). This shows that the point (t* = 0, s* = 0, m* = 
0) constitutes an infrared-stable fixed point with vanishing physical string 
tension T. This point is characterized by long-range correlations g{d) = 
27r^/a, with a constant a, and by the scaling law 

dl = ^p*d\ p*^(l-^) , (6.14) 

which shows that the Hausdorff dimension of world sheets is D^ = 2. For 
s = 0, the constant a can be computed analytically: 



2 ^. 1 + (2t/m2c2) + Vl + 4t/m2c2 
a = lim , (6.15) 

from which we recognize that 1 < p* < 2. 

At the infrared fixed point we can remove the cutoff. The renormaliza- 

tion of the model is easily obtained by noting that the effective action for 

transverse fluctuations to quadratic order decouples from other modes and 

is identical with the second term in (6.3) with D^ = d'^/p and p taking its 

saddle-point value. From here we identify the physical tension, stiffness and 

mass as 

T = A^{t + l), S = -, M = Amp. (6.16) 

P 

For s = 0, we can compute analytically the corresponding 7-functions: 

.-^lnt = 2 + of- 
a A \m 

,.„.-aA,„,„ = i + o(|j). (0.18) 

The vicinity of the infrared fixed point defines a new theory of smooth 
strings for which the range of the orientational correlations in embedding 



7t =-A— lnt = 2 + 0( — ) > (6-17) 
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space is always of the same order or bigger than the length scale 1/vT 
associated with the tension. The naively irrelevant term V^/M'^ in (6.1) 
becomes relevant in the large-D approximation since it generates a string 
tension proportional to M^ which takes over the control of the fluctuations 
after the orientational correlations die off. Note moreover, that it is exactly 
this new quartic term which guarantees that the spectrum p'^V (T, s, M, L, p) 
has no other pole than p = 0, contrary to the rigid string, which necessarily 
has a ghost pole sX p^ = —T/ s. 

By studying the finite-size scaling [87] of the Euclidean model (6.1) on a 
cylinder of (spatial) circumference R it is possible to determine the univer- 
sality class of confining strings. In the limit of large R, the effective action 
on the cylinder takes the form [41] 

S^^ ^ TTc(D-2) 

lim —=TR-^-—^ + ... , 6.19 

/3^oo p 6R 

for (D—2) transverse degrees of freedom, the universality class being encoded 
in the pure number c. This suggests that the effective theory describing the 
infrared behavior is a conformal field theory with central charge c. In this 
case the number c also fixes the Liischer term [88] in the quark-antiquark 
potential: 

n^)='ri?-^^^^^ + ... . (6.20) 

In [41] it has been shown that confining strings are characterized by c = 1. 
Although they share the same value of c, confining strings are clearly different 
(c = l)-theories than Nambu-Goto or rigid strings. Indeed, the former are 
smooth strings on any scale, while the latter crumple and fill the ambient 
space, at least in the infrared region. Our result c = 1 is in agreement with 
recent precision numerical determinations of Liischer and Weisz [30] of this 
constant. 

Having established that the model (6.1) describes smooth strings with 
c = 1, the question arises as to how much these results depend on the 
truncation of the original non-local action after the P^-term. In [41] it has 
been proved that the answer to this question is no: the value of c and the 
smooth geometric properties are independent of an infinite set of truncations, 
provided that a solution for the polynomial "gap" equation exists. These 
properties are presumably common to a large class of non-local world-sheet 
interactions. 
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7. High-temperature behavior of confining strings 

The high-temperature behavior of large- A^ QCD has been studied by Polchin- 
ski in [46], where he shows that the deconfining transition in QCD is due 
to the condensation of Wilson lines, and the partition function of QCD 
flux tubes can be continued above the deconfining transition; this high- 
temperature continuation can be evaluated perturbatively. So, any string 
theory that is equivalent to QCD must reproduce this behavior. However, 
the Nambu-Goto action has the wrong temperature dependence, while the 
rigid string parametrically has the correct high-temperature behavior but 
with the wrong sign and an imaginary part signaling a world-sheet instabil- 
ity [89]. 

The high-temperature behavior of the confining-string model proposed 
in [40] has been studied in [45], where it has been shown that this model has 
a high-temperature behavior that agrees in temperature dependence, sign 
and reality properties with the large- A^ QCD result [46]. The starting point 
will be again (6.2). In the Gauss' map, the value of the periodic coordinate 
^0 is — /3/2 < ^0 ^ /5/2 with (3 = 1/T and T the temperature. Note that, 
at high temperatures (/3 <C 1), the scale M^ can be temperature-dependent. 
This is not unusual in closed string theory as has been shown by Atick and 
Witten [90]. The value of ^i is -R/2 < ^i < R/2; ct>'{0's describe the 
D — 2 transverse fluctuations. We look for a saddle-point solution with a 
diagonal metric gab = diag (po ) Pi ) > arid a Lagrange multiplier of the form 
A"^ = diag(Ao/po,Ai/pi). 

After integration over transverse fluctuations we obtain, in the limit R -^ 
oo, four gap equations: 

^—^ = 0, (7.1) 

Po 

1 D-2 1 V2M , , 

px 2 4/3 (Ai + 1)3/4 ' ^ • ^ 



1/ , X 1 /, N , -D-2 vr 



(i-Ai) + ^^(Ai+t)-t-Ao + ^^:^ = 0, (7.3) 



1 /, ^ -D-2 1 

(t-Ai) (Ai+t) + 



vr 



Pi ' 2/3 

and a simplifled form of the effective action: 



2M(Ai+t)^/^ ^ 



0, (7.4) 



Z-eff 




^extTW- , (7.5) 
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with D = 2(Ai + t) representing the physical string tension. By inserting 
(7.2) into (7.4), we obtain an equation for (Ai + 1) alone: 

(A, + t) - ^1^ V2M (A, + *) VV ^^ - * = . (7.0) 

Without loss of generality we set 

(Ai + t)V^ = ^^, (7.7) 

7 

where 7 is a dimensionless parameter. It is possible to show that, at high 
temperatures, when 

tP' « ^^ , (7.8) 

we can completely neglect t in (7.6). Indeed, as we now show, Ai is propor- 
tional to (D — 2)^//3^. We can thus rewrite (7.6) as 



D-2 5 1/2 D-2 vr 



Ai 7;lK + ^ = 0. (7.9) 

2 8/3'^ o o«2 y ! 

We now restrict ourselves to the regime 



G{1, D) = f^l' (^) ' - 2. ^ > , (7.10) 



for which (7.9) admits two real solutions: 



5 D-2 1 



(A2)1/2 = ^ 7 -VG(7,Z)) . (7.12) 

V 1; 16/3 ' 2 2[3^ ^" ' ^ ' 

In both cases, Ai is proportional to (-D — 2)//3^, which justifies neglecting t 
in (7.6) and implies that the scale M^ must be chosen proportional to l//3^. 
Moreover, since the physical string tension is real we are guaranteed that 
M^ > 0, as required by the stability of our model. Any complex solutions 
for T> would have been incompatible with the stability of the truncation. 

Let us start by analyzing the first solution (7.11). By inserting (7.11) 
into (7.3), we obtain the following equation for pi : 

1 4 

1-- .._ ... ■ (7.13) 



PI 5 + . /25 - 4f^ 
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Owing to the condition (7.10), l/pi is positive and, since Af is real, the 
squared free energy is also positive: 



^^(«-|f-^(li^¥-i^/^^(^) ('-i77 



5+J25- 



1287r 
.2 -D-2 
2 



'(7.14) 

In this case the high-temperature behavior is the same as in QCD, but the 
sign is wrong, exactly as for the rigid string. There is, however, a crucial 
difference: (7.14) is real, while the squared free energy for the rigid string is 
imaginary, signaling an instability in the model. 

If we now look at the behavior of pi at low temperatures, below the 
deconfining transition [41], we see that 1/pi is positive. The deconfining 
transition is indeed determined by the vanishing of 1/pi at f3 = /3dee- In the 
case of (7.11) this means that 1/pi is positive below the Hagedorn transi- 
tion, touches zero at (3dec and remains positive above it. Exactly the same 
will happen also for F^, which is positive below /3dec ) touches zero at /3dec 
and remains positive above it. This solution thus describes an unphysical 
"mirror" of the low-temperature behavior of the confining string, without a 
real deconfining Hagedorn transition. For this reason we discard it. 

Let us now study the solution (7.12). Again, by inserting (7.12) into 
(7.3), we obtain 



1 



Pi 5 - , /25 



In this case, when 



-K /D-2 



1287r 
' 2 



-1/2 



(7.15) 



1/pi becomes negative. The condition (7.16) is consistent with (7.10) and 
will be taken to fix the values of the range of parameter 7 that enters (7.7). 
We will restrict to those that satisfy (7.16). Since po = 1 and Ai is real 
and proportional to l//9^, we obtain the following form of the squared free 
energy: 



(7.17) 



25-41^ 
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In the range defined by (7.16) this is negative. For this solution, thus, both 
1/pi and F'^ pass from positive values at low temperatures to negative values 
at high temperatures, exactly as one would expect for a string model under- 
going the Hagedorn transition at an intermediate temperature. In fact, this 
is also what happens in the rigid-string case, but there, above the Hagedorn 
transition, there is a second transition above which, at high temperature, 
Ai becomes large and essentially imaginary, giving a positive squared free 
energy. This second transition is absent in our model. 

Let us now compare the result (7.17) with the corresponding one for 
large- A^ QCD [46]: 



2gHm 
7r2/32 



F^(/3)qcd = -^^4^ , (7.18) 



where g^{l3) is the QCD coupling constant. First of all let us simplify our 
result by choosing large values of 7: 



1287r fD-2 
7> 



^1/2 



25 

In this case (7.17) reduces to 

F^(/3) = -- — ^^ (719) 

This corresponds exactly to the QCD result (7.18) with the identifications 

2 1 

g Of ^ > 

T 
N ocD-2 . 

The weak /3-dependence of the QCD coupling g'^{/3) can be accommodated 
in the parameter 7. Note that our result is valid at large values of 7, i.e. 
small values of g'^, as it should be for QCD at high temperatures [91]. Note 
also the interesting identification between the order of the gauge group and 
the number of transverse space-time dimensions. Moreover, since the sign of 
Ai does not change at high temperatures, the field x^ is not unstable. The 
opposite happens in the rigid-string case [89], where the change of sign of Ai 
gives rise to a world-sheet instability. 
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8. The influence of matter flelds on the deconflnement phase 
transition in the 3D GG model at flnite temperature 

8.1. Introduction 

The phase structure of the 3D GG model at finite temperature has for the 
first time been addressed in Ref. [51], where it has been shown that, in the 
absence of W-bosons, the weakly coupled monopole plasma undergoes the 
Berezinsky-Kosterlitz-Thouless (BKT) [82] phase transition into the molecu- 
lar phase at the temperature Tbkt = 5^/(2vr). Then, in Ref. [52], it has been 
shown that the true phase transition, which takes into account W-bosons, 
occurs at approximately half this temperature. Let us first discuss in some 
more detail the above-mentioned BKT phase transition, which takes place 
in the absence of W-bosons (i.e. in the continuum limit of the 3D lattice 
compact QED extended by the Higgs field). 

At finite temperature T = 1/(3, equations of motion of the fields x ^^'^ 
■0, entering the partition function (2.2), 

^mon = / VxDtpexpl - / cfxC[x,ij\gm,C] 

should be supplemented by the periodic boundary conditions in the temporal 
direction, with the period equal to f3. Because of that, the lines of magnetic 
field emitted by a monopole cannot cross the boundary of the one-period 
region. Consequently, at the distances larger than /? in the direction perpen- 
dicular to the temporal one, magnetic-field lines approaching the boundary 
should run almost parallel to it. Therefore, monopoles separated by such 
distances interact via the 2D Coulomb potential, rather than the 3D one. 
Since the average distance between monopoles in the plasma is of the or- 
der of C~^ ) we see that, at T > C^ > the monopole ensemble becomes 
two-dimensional. Owing to the fact that C is exponentially small in the 
weak-coupling regime under discussion, the idea of dimensional reduction is 
perfectly applicable at the temperatures of the order of the above-mentioned 
critical temperature Tbkt- 

Note that, due to the T-dependence of the strength of the monopole- 
antimonopole interaction, which is a consequence of the dimensional reduc- 
tion, the BKT phase transition in the 3D Georgi-Glashow model is inverse 
with respect to the standard one of the 2D XY model. Namely, monopoles 
exist in the plasma phase at the temperatures below Tbkt and in the molec- 
ular phase above this temperature. As has already been discussed, the anal- 
ogy with the 2D XY-model established in Ref. [62] is that spin waves of that 
model correspond to free photons of the 3D GG model at zero temperature. 
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while vortices correspond to magnetic monopoles. At finite temperature, 
disorder is rather produced by the thermally excited W-bosons [52], whereas 
monopoles order the system, binding W-bosons into pairs. However, the 
analogy is still true in the case of the continuum version of 3D compact 
QED (with the Higgs field) under discussion. 

Let us therefore briefly discuss the BKT phase transition, occurring in 
the 2D XY-model at a certain critical temperature T = ?bkt. At T < 
-^BKT' ^^^ spectrum of the model is dominated by massless spin waves, and 
the periodicity of the angular variable is unimportant in this phase. The 
spin waves are unable to disorder the spin-spin correlation functions, and 
these decrease at large distances by some power law. On the contrary, at 
T > Tgj|rp, the periodicity of the angular variable becomes important. This 
leads to the appearance of topological singularities (vortices) of the angular 
variable, which, contrary to spin waves, have nonvanishing winding numbers. 
Such vortices condense and disorder the spin-spin correlation functions, so 
that those start decreasing exponentially with the distance. Thus, the nature 
of the BKT phase transition is the condensation of vortices at T > T^^rj,. 
In another words, at T > Tgj^rp, free vortices do exist and mix in the ground 
state (vortex condensate) of an indefinite global vorticity. Contrary to that, 
at r < Tgj^rp, free vortices cannot exist, and they rather mutually couple into 
bound states of vort ex- ant i vortex pairs. Such vortex-antivortex molecules 
are small-sized and short-living (virtual) objects. Their dipole-type fields 
are short-ranged and therefore cannot disorder significantly the spin-spin 
correlation functions. However, when the temperature starts rising, the size 
of these molecules grows, until at T = T^^r^ it diverges, that corresponds to 
the dissociation of the molecules into pairs. Therefore, one of the methods to 
determine the critical temperature of the BKT phase transition is to evaluate 
the mean squared separation in the molecule and to find the temperature at 
which it starts diverging. 

Let us now return to the continuum limit of the finite-temperature 3D 
compact QED, extended by the Higgs field, and determine there the mean 
squared separation in the monopole-antimonopole molecule. One can then 
see that, up to exponentially small corrections, the respective critical tem- 
perature is unaffected by the finiteness of the Higgs-boson mass. Indeed, the 
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mean squared separation reads' 



J a x|x| 9 exp 



L 



2\ _ l^l>™iv 



47rT 



Ko("^//|x| 



J (i^x|x| ~s^ exp ^2^Ko (^tt-hIxI) 



|x|>mv 



l.l) 



where Kq denotes the modified Bessel function. Disregarding the exponential 
factors in the numerator and denominator of this equation, we obtain (-^^) — 

i^io~'T — 2T) that yields the above-mentioned value of the BKT critical 
temperature Tbkt = 5^/(2vr). Besides that, we see that, as long as T does 
not tend to TbkTj and in the weak-coupling regime under study, the value 
of ■\/{L'^) is exponentially smaller than the characteristic distance in the 
monopole plasma, C~^ > i-6- molecules are very small-sized with respect to 
that distance. 

The factor (3 at the action of the dimensionally-reduced theory, S'd.-r. = 
(3 J CpxC[x-,4'\grniC,]i can be removed [and this action can be cast to the 
original form of Eq. (2.2) with the substitution d?x -^ d^x] by the obvious 
rescaling: 



5d.. 



d xC 



r''^,r 



K,f3C 



U) 



Here, K = gl^T, x""^"^ = VPx, V^ = \/?^, and in what follows we will 
denote for brevity x^^^^ and tp'^'^^ simply as x and tp, respectively. Averaging 
then over the field ip with the use of the cumulant expansion we arrive at 
the following action: 



-2e 



Sd.-r. 



d xd ycos 



d^x 



lid.X? 



2,^ cos 
!^x(x))/C(2)(x-y)cos 



Kx 



Kx{y) 



(8.3) 



In this expression, similarly to Eq. (2.6), we have disregarded all the cu- 
mulants higher than the quadratic one, and the limits of applicability of 
this bilocal approximation will be discussed below. Further, in Eq. (8.3), 



x;(2)(x) 



^KDi^li^) 



1, where I^^U^) = i^o(m-i?|x|)/(27r) is the 2D 



Yukawa propagator, and ^ = P( exp 



K r^{2) 
2 ^mn 



(W) 



denotes the monopole 



fugacity modified by the interaction of monopoles via the Higgs field. The 



' In this section, 3- vectors are denoted as a, whereas 2- vectors are denoted as a . 
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exponential factor entering ^ reads 



^ oc exp 



4vr / _, /eT 

— ( m-iye + Tin ( — c 



Here, we have introduced the notation c = niH/rnw, c < 1, and 7 ~ 0.577 is 
the Euler constant, so that ^ ~ 0.89 < 1. We see that the modified fugacity 
remains exponentially small, provided that 



T < 



mw^ 



In 



fc-' 



(8.4) 



This constraint should be updated by another one, which would pro- 
vide the convergence of the cumulant expansion. Analogously to the zero- 
temperature case, the divergence of the cumulant expansion would indicate 
that the Higgs vacuum loses its normal stochastic property and becomes a 
coherent one. In order to get the new constraint, notice that the parameter of 
the cumulant expansion reads ^I^'^', where I^'^' = J (Px)C^'^\x.). Evaluation 
of the integral /^^^ yields [53] : 



/(2) ^ 



2tt 



m 



H 




) 87rT 
~^ 1- 



SttT 



4ttT 

9'^ 



+ e^ 



1 + 



(8.5) 



(Note that, at T 



91 Iz 

47r' 



SttT 
9^ 



4ttT 



-21nc, i.e. /'^' remains finite.) In the 



derivation of this expression, the parameter a = A'K\f2/nT j g^ was assumed 
to be of the order of unity. That is because the temperatures we are working 
at are of the order of Tbkt- Due to the exponential term in Eq. (8.5), the 
violation of the cumulant expansion may occur at high enough temperatures 
[that parallels the above-obtained constraint (8.4)]. The most essential, ex- 
ponential, part of the parameter of the cumulant expansion thus reads 



^/^ ' ex exp 



47r 



m-we -I- T 



In 



Therefore, the cumulant expansion converges at the temperatures obeying 
the inequality 



T < 



2-K 



In (f c) 



which strengthens the inequality (8.4). On the other hand, since we are 
working in the plasma phase, i.e. T < TbkT) it is enough to impose the fol- 



February 1, 2008 18:47 WSPC/Trim Size: 9.75in x 6. Sin for Proceedings Antonov 



54 D. Antonov and M. C. Diamantini 

lowing upper bound on the parameter of the weak-coupling approximation, 
k: 

27re 



K< 



'^ -ln(fc) 



Note that, although this inequality is satisfied automatically at ^c ~ 1 (or 



c ~ 1), since it then takes the form k < v27ree, this is not so automatic in 
the BPS limit, c <C 1. Indeed, in this case, we have kIu (^) < 27re. It is, 
however, quite feasible to obey this inequality, since the logarithm is a weak 
function. 

8.2. Higgs-inspired corrections to the RG flow in the absence 
of W-bosons 

Although we have seen in the previous subsection that the propagating Higgs 
field does not change the value of the critical temperature Tbkt, it is instruc- 
tive to derive the corrections it produces to the RG flow. Such a derivation 
can also be extended to the SU(A^)-case (2.5). In that case [83], it is quali- 
tatively clear that Tbkt should remain the same, since it only diff'ers from 
that of the SU(2)-case by the factor ^^, which is equal to unity. However, 
some peculiarities of the SU(A^)-case at A^ > 2 become clear only upon a 
derivation of the RG equations, which will be given below. 

Let us thus start with the SU(2)-case (2.2), whose action after the di- 
mensional reduction has the form (8.2). In what follows, we will use the 
usual RG strategy based on the integration over the high-frequency modes. 
Note that this procedure will be applied to all the flelds, i.e. not only to 
X, but also to ip. Splitting the momenta into two ranges, < p < A' and 
A' < p < A, one can define the high-frequency modes as h = XA — XA', 

4> = i^A — V'A') where Oa'(^) = / 72^^*^^^(P) ^^^ consequently, e.g. 

0<p<A' 

/i(x) = J ,„ ?i e^P^x(p). The partition function, 
A'<p<A 

^A = / Vxip)Vtp{p) exp {-S'd.-r. [XA, V'a]} , 
0<p<A 
can be rewritten as follows: 

Za = y" Vx{p)ViP{p) exp 1^ y" (fx [XA'S^XA' + ^A' {d^ - m%) V'a'] | Z\ 
0<p<A' 
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where 



A'<p<A 



Px(p)W'(p)e.p{//. 

+2ee^^(^A'+<^) cos (^ (XA' + h)M\ 



\hd^h + ]^<p {d^ - ml) 0+ 



and ^ = f3C,. Owing to the exponential smallness of the fugacity, Z' can 
further be expanded as 

Z'~l + 2e /'(i2x(e^(^A'+<^)^ ^cos(/^(XA' + /i))) +2f fd^xd^yx 

,v^(V'A,(x)+<^(x))gVi?(^A'{y)+<^(y))^ ^cos (Vk (xA'(x) + /i(x))) X 
xcos(/^(XA'(y) + /i(y)))) -(e^(^A'W+^W)) ^e^('^A'{y)+<^(y))^ x 



X ( cos ( VK (xA'(x) + /i(x)) ) )^ ^cos (^Vi^ (xA'(y) + Hy))^ I ^ 



where 



kO), 



J Vx{p)exp{^jd^xhd^h)0 

A'<p<A 

f Vx{p)exp{^fd'^xhd^h) 
A'<p<A 



'^ m? 



O 



(O), 



/ DV'(p)exp [i Jd^x0 (3' — iii'H 
A'<p<A 

f PV(P) exp [i / d'^x^ (92 - m|^) (/)] 
A'<p<A 



Carrying out the averages we arrive at the following expression for Z': 
Z'^1^ 2,AmiO) I d^xe^^^' cos ( V^x.) + (^^(0)^.(0))^ / d^xd 
^gV^(^A,(x)+VA.(y)) ^ J^2fc(^ _ y)^2(^ - y) - l" 

k=±l 

X cos [\/^(xa'(x) + A:xA'(y)) 
where A(x) = e-^^^W/^ 5(x) e= e^^^W/^, 



t/x 



(8.6) 



GhU) 



d'^p 



a«PX 



A'<p<A 



(27r)2 p 



2 n2 



G^fx) 



d^p e' 



ipx 



A'<p<A 



(27r)2 p2 _|_ jyij^ 
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Since in what follows we will take A' = A — dA, the factors 
[^^ (x — y)i?^(x — y) — l] , /c = ±1, are small. Owing to this fact, it is con- 
venient to introduce the coordinates r = x — y and R = 2 (x + y ) , and Taylor 
expand Eq. (8.6) in powers of r. Clearly, this expansion should be performed 

up to the induced-interaction term ~ ^^^ / (PRe^^^"^^' ^^^ cos ( 2\/Kxa' (R-) ) ; 
that itself should already be disregarded. As a result, we obtain the following 
expression for Zj\^ (For the sake of uniformity, we replace d?R by cPx): 

Za= y Px(p)^V'(p)exp|/'d2a;UvA'(5^-"^H)^A' + 



0<p<A' 



+a2{iA{Q)B{^)Ye 



2 2v^Va' _ i 



1 + ai(e^(0)5(0))2:^e2^'^A' 



+2C^(0)B(0)e^'^A' cos (/^XA') }} , 
where we have introduced the notations 

ai = f d^rr^ [A-\r)B^{r) - l] , 02 = / dV [A-\r)B\r) - 1 

Taking into account that A' = A — dA it is straightforward to get 



{d^XA'f + 
(8.7) 



dA f A2 
ai = aiK—^ ( 1 H j- 



A5 



m 



H 



dA / A2 
02 = a2K-p^ [ 1 H 2" 



A3 



m 



H 



Here, ai^2 stand for some momentum-space-slicing dependent positive con- 
stants, whose concrete values will turn out to be unimportant for the final 
expressions describing the RG flow. 

Next, since 01^2 occur to be infinitesimal (being proportional to dA), the 
terms containing these constants on the r.h.s. of Eq. (8.7) can be treated 
in the leading-order approximation of the cumulant expansion that we will 
apply for the average over ip. In fact, we have 



Vil){p) exp 



d^xil^^i {p^ — rn^j) Va' 



0<p<A' 



exp (b f d^xe^^^^'fix)] ~ exp fte^^^W f d^xfix) + 



H e 

2 



4E'G{0) 



y"d2^A(e^^«(-y)-l)/(x)/(y) 
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where / is equal either to unity or to {d^XA')'^, b ~ ai^2{^A{0)B{0))'^, and 

cPp 



G(x)^ 



=«PX 



0<p<A' 



(27r)2p2 + ^2^' ^(0) 



In 1 + 



in 



niH^ 



In order to estimate the parameter of the cumulant expansion, k = 
^g2_ft'G(0) J ^2^ ^g4_ft'G(x) _ 2^ ^ note that we are working in the phase where 
monopoles form the plasma, i.e. below Tbkt- Because of this fact, 
4i^|G(x)| < 32TT{A'/mH)'^, that, due to the factor {A'/mn)'^, is generally 
much smaller than unity. Owing to that, we get 



K ~ 6 1 + 



A' 



/2 



rriH" 



■AK f(fxG{x) ~ 



2bK 



Trm 



H 



1 + 



K ( h! 



2tt \mH 



Choosing for concreteness b = 02(^^(0)5(0))^ and taking into account that 



Am . 1 - -0,(0) = 1 - j^x • 

K ^ , , K f A V dk 

i^(0).l + -G,(0).l + -(— ) - 



(8.9) 
(8.10) 



we obtain to the leading order: k ~ 2a2K ^"^ / {-Km?^) . This quantity possesses 
double smallness - firstly, because 02 is infinitesimal and, secondly, due to 
the exponential smallness of ^. 

Such an extremely rapid convergence of the cumulant expansion enables 
one to replace q^^-^'Va' in the terms proportional to ai^2 on the r.h.s. of 

. Com- 



Eq. (8.7) by the average value of this exponent equal to I 1 H j 

paring the so-obtained expression with the initial one, we arrive at the fol- 
lowing renormalizations of fields and parameters of the Lagrangian: 



new /^ / new /^ / j^n 

Xa' =Cxa', ^a' =C'V'a', K 



K 

C2' 



/^ 



Q2 



, e 



^(o)i?(o)e, 

(8.11) 



where /U 



m 



_ff' 



C 



1 + :^(^^(o)s(o))2 

Kax 



1 + 



1 + 



(^^(0)S(0))2 1 + 





1/2 


A2 K\-] 


1/2 


m?^ 2n ) _ 





^.12) 



February 1, 2008 18:47 WSPC/Trim Size: 9.75in x 6. Sin for Proceedings 



Antonov 



58 



D. Antonov and M. C. Diamantini 



Besides that, we obtain the following shift of the free-energy density F 

\nZ' . 



F 



a2(e^(0)B(0))2 1 + 



A 



K 
/2 \ 27r 



rriH^ 



F" 



a2(^A(0)5(0))Ml + 



K A" 
2-K mu"^ 



^.13) 



where V is the 2D-volunie (i.e. area) of the system. 

By making use of the relations (8.9), (8.10), it is further straightforward 
to derive from Eqs. (8.11)-(8.13) the RG equations in differential form. These 
are 



^-^^-f)f. — f-'^^ 



An 



1+ — + 1 

27r 



A2 



dk 
A5 



d/x = -^(KOV 






dh 

A5 



dF = QsKr 






dh 

A3 



Now let us change from the momentum scale to the real-space one: A 
a = 1/A, dK -^ —dA. This modifies the above equations to 






4-K a 



dK 



dfi 



-^K^ea^da 



-^(KO'f^a'da 



K 

1+ hi 

V2vr 

fK 

l+\ — + 1 
V2vr 



/ia^ 



/iO^ 



dF = a2K^ ada 



1 + 



K 

2^ 



+ 1 



1 



yua^ 



^.14) 
^.15) 
^.16) 

^.17) 



Our main aim below is to derive from Eqs. (8.14)-(8.16) the leading- 
order corrections in {^a?)^^ to the BKT RG flow in the vicinity of the 



'(0) 



(0) 



critical point, Kc = Svr (that clearly corresponds to Tbkt); ^c = 0, where 
y = ^o?, and the subscript "'''' " denotes the zeroth order in the {na?')~^- 
expansion. These values of Kc and yc will be derived below. Additionally, 
it will be demonstrated that, in the critical region, ^ is evolving very slowly. 
Owing to this fact, the initial assumption on the largeness of ^ (namely, 
that it is of the order of m^^), will be preserved by the RG flow, at least in 
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that region. This enables one to consider fj, as almost a constant and seek 
corrections to the RG flow of K^^' in powers of (fia^)^^. The zeroth-order 
equation stemming from Eq. (8.15) then reads 

dK^^^ = -^K^^^^ea^da. (8.18) 

Respectively, the zeroth-order in {pio?)^^ equation for y has the form 

dy(0)2 = 2— y^V (8.19) 

a 

where a; = 2 ^^. Equations (8.18) and (8.19) yield the above-mentioned 

leading critical value of K, Kc ■ Next, with this value of Kc , Eq. (8.18) 
can be rewritten in the vicinity of the critical point as 

dx = {STrfaifa^da. (8.20) 

Introducing further instead of y the new variable z = (87r)^aiy^ and per- 
forming the rescaling a^'^'" = ay87rai/a2 we get from Eqs. (8.17), (8.19), 
and (8.20) the following system of equations: 

dzio) = 2-xz(0), dx = z(o)^, dF(0) = ^(0)^. (g^l) 

a a a-^ 

These equations yield the standard RG flow in the vicinity of the critical 
point, X = z^^' = 0, which has the form [82, 84] z^^> — x^ = t, where 
r oc (Tbkt — r)/2BKT is some constant. In particular, x ~ ^/z at T = 

/ /g-|\— 1/2 

Tbkt — 0. Owing to the first of Eqs. (8.21), this relation yields ( z)^J J — 

1/9 

(z^'^') = ln(a/ain), where the subscript "m" means the initial value. 

Taking into account that z\J is exponentially small, while z^^' ~ 1 (the 
value at which the growth of z^^' stops), we obtain in the case xi^ ^ \/^- 

Zjjj J ~ r~^". According to this relation, at T = Tbkt — 

0, the correlation length diverges with an essential singularity as aij^ ~ 
exp (const /-yr). (At T < TbkT) a = d, while at T > TbkT) the correlation 
length becomes infinite due to the short-rangeness of molecular fields.) As 
far as the leading part of the free-energy density is concerned, it scales as 
i^(0) rs^ Q-2 g^j^j therefore remains continuous in the critical region. Moreover, 
the correction to this behavior stemming from the finiteness of the Higgs- 
boson mass [the last term on the r.h.s. of Eq. (8.17)] is clearly of the same 
functional form, ~ exp (— const'/\/r), i.e. it is also continuous. 

We are now in the position to address the leading-order [in (fia^)^^] cor- 
rections to the above-discussed BKT RG flow of K^^> and z^^' . To this end. 
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let US represent K and z as K = K^^'^+K^'^'^ /{jia?), z = z^^'^ + z^'^'^ / {jxa?) that, 
by virtue of Eqs. (8.14) and (8.15), leads to the following novel equations: 

.,.<« _ 2A-<«^ = -4,^ (.<«. + £!^ Vl + ^ + f:;^) . (8,22) 

<i.-(l)-2.'"^ = -2^ 2wf:^-2l+ — Ca-M-A-'"') , (8,23) 

In the vicinity of the critical point, we can insert into Eq. (8.23) the above- 
obtained critical values of K^^' and z^^' to get 

dz(i)=2zW— . (8.24) 

a 

Therefore, z^^> = C\c? ^ where C\ is the integration constant of dimensional- 
ity (mass)^, C\ <^ fj,. Inserting further this solution into Eq. (8.22), consid- 
ered in the vicinity of the critical point, one obtains the following equation: 

<,A-<')-2A-<')^ = -l^^('?^ + 5V (8,25) 

a fi a \ Sir I 

Its integration is straightforward and yields 

Aa..C,(,.»V-*.J2^, (8,2.) 

where the dimensionless integration constant C2 should be much smaller 

than (fJ-a^) "''' . [Note that the last term in Eq. (8.26) is positive.] There- 
fore, the total correction, K^^> /{fia^), approximately scales with a in the 
critical region as (4 _^c ) a'^ ' ^® ®^^ that, at the critical point, this ex- 
pression vanishes due to the divergence of the correlation length. Note also 
that, for z'^^ <^ /ia^, Eq. (8.26) can obviously be rewritten as the following 
dependence of K^^' / (fia^) on z'^': 

With the above-discussed critical behavior of the correlation length, a^r), 
this relation determines the correction to the BKT RG flow, z^^' — 

2-/e)'=- 

Finally, in order to justify the approximation that /x was considered con- 
stant, we should check that, under the RG flow, ^ indeed evolves slowly. To 
this end, let us pass in Eq. (8.16), considered in the critical region, from the 



February 1, 2008 18:47 WSPC/Trim Size: 9.75in x 6. Sin for Proceedings 



Antonov 



3D Georgi-Glashow model and confining strings 



61 



variable ^ to the above-introduced variable z and again perform the rescal- 
a"^™. This yields ^ = -f^ord^ = -^|i. Since Ci < /i, 



mg a 



M 



da 



we conclude that -L^ <^ ^. The obtained inequality means that, in the 
vicinity of the BKT critical point, /x is really evolving slowly. This fact jus- 
tifies our treatment of /x as a large (with respect to A^) constant quantity, 
approximately equal to its initial value, of the order of m|^. 

Let us now proceed with a similar RG analysis of the SU(A^)-version of the 
theory (2.2) given by Eq. (2.5). Applying to the respective dimensionally- 
reduced theory the above-described RG procedure, we arrive at the following 
analogue of Eq. (8.7): 



Za= y Vxip)V^P{p)explJ(fxl^i;A>{d^-ml)iPA,+ 

K{q^-qJ)xA'] +2^A(0)5(0)x 



0<p<A' 



Oa COS 



u 



'^e^^A' J2 COS (VKqa A 



5"^ + (e^(0)S(0))2— e^^^^A' : 
8 



XI "i' (^« + ?jT (.Qi + Qj)^ cos VK {qi - qj) xa 



ij 



d,.XA'j[d,.XA' 




Here, a, (3 = 1, . . . , (A^ — 1), and we have introduced notations similar to 



a-^ = Id rr 



B'^M^Kq^qjChir) _ -^ 



, a2 = I d'^r 



so that, at A' = A - dA, 



^j 



A2 



A*^ V -^ mj 



B'^M^Kq^qjChir) _ ^ 



A2 



, U2 



r.dA 



The main difference of Eq. (8.27) from Eq. (8.7) is due to the terms contain- 



ing COS 



K[ 



It - Qj) XA' 



which violate the RG invariance. Nevertheless, 
approximately this invariance does hold, since the respective sums are dom- 
inated by the terms with i = j. Working within this approximation and 
making use of the identity ^ qfQi = "j"*^ ' ^^ obtain 



Za ~ exp 



02 



m^^^^Ammf U^ '' 



K 
2 \ 2^ 



ruH 



V 



February 1, 2008 18:47 WSPC/Trim Size: 9.75in x 6. Sin for Proceedings 



Antonov 



62 



D. Antonov and M. C. Diamantini 



J Vx{p)V%b{p) exp \ Jd^xl ^V'A' {d^ - m]j) ^a'+ 



0<p<A' 



+2^y4(0)5(0)e^-^'^A' ^cos i^qiXk 



1+ 1 + 



K 
^,2 \ 2. ^^^^ 



rriH^ 



ia{0)B{0)f 



(^mXa 



- \2 



^.28) 



This expression has again been derived in the leading order of the cumulant 
expansion applied in course of the average over -0. 

The shift of the free-energy density stemming from Eq. (8.28) [cf. 
Eq. (8.13)] reads 



F" 



N(N -1) , ,, , . ^^9 / K A'^ 



27r rriH^ 



^.29) 



As far as the renormalization of fields and coupling constants is concerned, 
it is given by Eq. (8.11), where the first equation should be modified as 
^^^new _ Cx\/, and the parameter C from Eq. (8.12) now reads 



C 



i + ^mo)B{o)f 



i + ^mo)B{o)f 



1 + 



A' 



K-i 1/2 

2X2;;^ 



mn" 



1 + 



A2 K 



mjj 2iT 



1/2 



^.30) 



Prom Eqs. (8.29) and (8.30) we deduce that, in the SU(Af)-case at A^ > 2, the 
RG flow of couplings and of the free-energy density is identical to that of the 
SU(2)-case. Indeed, all the A'-dependence can be absorbed into the constants 
Q^i,2 by rescaling them as oti = Nai/2, a^ = N{N — 1)02/2 and further 
redefining [cf. the notations introduced after Eq. (8.20)] z = (87r)^aiy^ and 
^new _ a y/gT]- 0-^/(52. In particular, the critical temperature Tbkt remains 
the same as in the SU(2)-case. (As has already been mentioned, this also 
follows from the estimate of the mean squared separation in the monopole- 
antimonopole molecule, if one takes into account that the square of any 
root vector is equal to unity.) Thus, the principal difference of the SU(A^)- 
case, N > 2, from the SU(2)-one is that, while in the SU(2)-case the RG 
invariance is exact (modulo the negligibly small higher-order terms of the 
cumulant expansion applied to the average over ip), in the SU(A^)-case it is 
only approximate, even in the limit rriH 



00. 
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8.3. Finite-temperature 3D compact QED with massless 
fundamental fermions 

Let us consider the extension of the model (2.1) by fundamental dynami- 
cal quarks [55] (for the sake of simplicity, we omit the summation over the 
flavor indices, although consider the general case of an arbitrary number 
of flavors). The action is modified by AS" = —i J d^x-ip i^D + /i^$" 1 xp, 

where the Yukawa coupling h has the dimensionality [mass]^'^, Df^ip = 
{d^j_ — ig^A'f^j Tp, tl} = Tp'P, the Euclidean Dirac matrices are defined as 

7 = —ij3a with /3=(„ ),a=( w„)7 ^'^^ '^ denote the Pauli matrices. 

Vo-i; \tOJ 

As will be demonstrated, at T > TbkT) quark zero modes in the monopole 
field lead to the additional attraction of a monopole and an antimonopole 
(M and M for shortness) in the molecule. In particular, when the number of 
these modes (equal to the number of massless flavors) is sufficiently large, the 
molecule shrinks so strongly that its size becomes of the order of the inverse 
W-boson mass. Another factor which defines the size of the molecule is the 
characteristic range of localization of zero modes. Namely, it can be shown 
that the stronger zero modes are localized in the vicinity of the monopole 
center, the smaller the molecular size is. We will consider the case when the 
Yukawa coupling h of quarks with the Higgs field vanishes, and originally 
massless quarks do not acquire any mass. This means that zero modes are 
maximally delocalized. Such a weakness of the quark-mediated interaction 
of monopoles opens the possibility for molecules to undergo eventually the 
phase transition into the plasma phase at the temperatures of the order of 
Tbkt- However, this will be shown to occur only provided that the num- 
ber of flavors is equal to one, whereas at any larger number of flavors, the 
critical temperature becomes exponentially small. This means that the in- 
teraction mediated by such a number of zero modes is already strong enough 
to maintain the molecular phase at any larger temperature. 

For a short while, let us consider the general case h ^ 0. One can then 
see that the Dirac equation in the field of the third isotopic component 
of the 't Hooft-Polyakov monopole [27] decomposes into two equations for 
the components of the SU(2)-doublet ip. The masses of these components 
stemming from such equations are equal to each other and read mq = hrj/2. 
Next, the Dirac equation in the full monopole potential has been shown [85] 
to possess a zero mode, whose asymptotic behavior at r = |x| ^ "^(7^ has 
the form 

Q-mqT 
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Here, Xn ^^^ the upper and the lower components of the mode, i.e. ip = (^-) . 
Next, n = 1,2 is the isotopic index, z^ = 1,2 is the Dirac index, s"*" = (g), 
s~ = (^), and M is the normahzation constant. 

It is a well known fact that, in 3D, the 't Hooft-Polyakov monopole is 
actually an instanton [3,26]. Therefore, we can use the results of Ref. [86] on 
the quark contribution to the effective action of the instanton-antiinstanton 
molecule in QCD. Let us thus recapitulate the analysis of Ref. [86] adopting 
it to our model. To this end, we fix the gauge <l>" = r]6^^ and define the 
analogue of the free propagator 5*0 by the relation Sq = —i i^d + niqT^ j . 
Next, we define the propagator Sm in the field of a monopole located at the 
origin, 1"^ [Af^ -^ e^^^xV {gr^) at r > m^^], by the formula S^^ = Sq^- 
g^^A"-^ . Obviously, the propagator S*^ in the field of an antimonopole 
located at a certain point R, A""^ (x) = —A""^ ix — R], is defined by the 



equation for S^ with the replacement A"" -^ A"" . Finally, one can 
consider the molecule made out of these monopole and antimonopole and 
define the total propagator S in the field of such a molecule. A"" = A""^ + 
A"-M ^ by means of the equation for S^ with A"-^ replaced by A'^. 

One can further introduce the notation |V'n), n = 0,1,2, .. ., for the eigen- 
functions of the operator —i^D defined at the field of the molecule, namely 
—i^D \tpn) = A„ \ipn), where Aq = 0. This yields the following formal spectral 
representation for the total propagator S: 



oo 



S {x, y) 



Sr^ \i^n{x)) (^n(^)| 

^ A„ - imqT^ 



Next, it is worth recalling the mean- field approximation, discussed in sec- 
tion 3, within which zero modes dominate in the quark propagator, i.e. 

Si^,^^)^\Mmi^^S,ix,y). (8.32) 

The approximation (8.32) remains valid for the molecular phase near the 

phase transition (i.e. when the temperature approaches Tbkt from above), 

we will be interested in. That is merely due to the fact that, in this regime, 

molecules become very much inflated being about to dissociate. 

/ _-, _-, _-, \ -1 

Within the notations adopted, one now has S = I Sj^ + S^, — Sq 
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SffS ^Sm, where 



;M\/ ;M \ih^^\/ih^\ 

O — Dq — (Dm — 'JQ) -Jq WM ~ '-'Oj — '-'0 : 3~ '-'o 



and IV'o )) IV'o ) ^'^^ ^^^ zevo modes of the operator —i^D defined in 
the field of a monopole and an antimonopole, respectively. Denoting fur- 
ther a = (V'o^|fl'7^^''^|V'o^)i if is straightforward to see, by the defi- 
nition of the zero mode, that a = (V'o^l - i^ d]^)^^) = (^o^|5'o"^|V'o^). 
This yields S = Sq + (a*/m^)|V'o^)(V'o^|) where the star stands for com- 
plex conjugation, and therefore det5 = [l + {\a\/mq) ] • detS'o- Finally, 
defining the desired effective action as F = ln[det 5~^/det S'jj' ], we obtain 
F = const + Nf In {m^ + jap) in the general case with Nf fiavors. The con- 
stant in this formula, standing for the sum of effective actions defined at the 
monopole and at the antimonopole, cancels out in the normalized expression 
for the mean squared separation in the MM-molecule. 

Let us now set h equal to zero, and so niq is equal to zero as well. No- 
tice first of all that, although in this case the direct Yukawa interaction 
of the Higgs boson with quarks is absent, they keep interacting with each 
other via the gauge field. Owing to this fact, the problem of finding a 
quark zero mode in the monopole field is still legitimate. [Note that, ac- 
cording to Eq. (8.31), this mode will be non-normalizable in the sense of 
a discrete spectrum. However, it is clear that, in the gapless case niq = 
under discussion, the zero mode, which lies exactly on the border of the 
two contiguous Dirac seas, should be treated not as an isolated state of 
a discrete spectrum, but rather as a state of a continuum spectrum. This 

means that it should be understood as follows: |^o(^)) ~ li^i (e*P''/r), where 

p— >o 

p = \'p\. Once considered in this way, zero modes are normalizable by the 
standard condition of normalization of the radial parts of spherical waves, 

oo 

Rpi, which reads [92] f drr^RpuRpi = 2it5{p' — p)] The dependence of 



the absolute value of the matrix element a on the distance R between a 
monopole and an antimonopole can now be readily found. Indeed, we have 
|a| oc J Sr / [r'^\r — n\\ ~ — 47rln(;ui?), where /U stands for the infrared cutoff. 
We can now switch on the temperature are explore a possible modifica- 
tion of the BKT critical temperature, TbkT) due to the zero-mode medi- 
ated interaction. As has been discussed, this can be done upon the eval- 
uation of the mean squared separation in the MM-molecule and further 
finding the temperature below which it starts diverging. In this way we 
should take into account that, in the dimensionally-reduced theory, the 
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usual Coulomb interaction of monopoles (Without loss of generality, we 
consider the molecule with the temporal component of R equal to zero.) 

1/2 

goes over into —2Tln{fiTZ), where IZ denotes 



+ 00 

R-^= J2 {n^ + {i3nf 

n=—oo 



the absolute value of the 2-vector R. This statement can be checked e.g. by 
virtue of the Euler - MacLaurin formula. As far as the novel logarithmic 



interaction, \n{pR) 
forms into 



+00 



Y: In fi{n' + {l3n 



,2\l/2' 



, is concerned, it trans- 



ttTTZ + In [1 - exp(-27^^7^)] - ln2 . 



.33) 



Let us prove this statement. To this end, we make use of the following 
formula [93]: 



E 

n=l 



1 



n^ + x^ 



1 

2^ 



7rcoth(7ra:) 



This yields 



+ 00 

^E 

n=—oo 



1 



-I 2 °o 

1 xa Y^ 



1 



a fax\ 

r, — - coth — . 

2 2 V 2 y 



x2 + {2Trn/af x 27r2 ^ n2 + (xa/27r) 
On the other hand, the l.h.s. of this expression can be written as 



1 a Y^ 

2dx ^ 

n=— oo 



In 



x' + 



27rn 



Integrating over x, we get 
In 



+ 00 



n=— oo 



X^ + 



27rn 



a / dxcoth 



ax 



2 In sinh 



ax 



ax + 21n(l-e"""=) -21n2 . 



Setting 



2n 



fj,f3 and x = fiTZ we arrive at Eq. (8.33). 



Thus, the statistical weight of the quark-mediated MM-interaction 
in the molecule at high temperatures has the form exp(— 2iVj In |a|) oc 
[ttTT?. + In [1 — exp(— 27rT7^)] — ln2]~ ^. Accounting for both, (former) 
logarithmic and Coulomb, interactions we eventually arrive at the follow- 
ing expression for the mean squared separation (L"^) in the molecule as a 
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function of T, g, and Nf. 

/ dTZTZ ~^ [ttTTZ + In [1 - exp(-27rT7^)] - In 2]"^^^ 
{L') = '^ —-, . (8.34) 



1"^ r ^-n , 1 n / o ^-nM 1 m-2N 



f dTZTZ ~^ [irTTZ + In [1 - exp(-27rT7^)] - In 2] 



/ 



'•w 



At large TZ, In 2 <C ttTTZ and |ln [1 - exp(-27^^7^)]| ~ exp(-27^^7^) < 
ttTTZ. Consequently, we see that {L"^) is finite at T > T^^rj, = (2 — 
Nf)g'^ /{4tt), that reproduces Tbkt at Nf = 0. For A'^^ = 1, the plasma phase 
is still present at T < (7^/(47r). Instead, when Nf > 2, at whatever tem- 
peratures, which are parametrically larger than the temperature of dimen- 
sional reduction, 0{C^'^), the monopole ensemble exists only in the molecu- 
lar phase. At such temperatures, massless dynamical quarks produce their 
own deconfincment. They also destroy the monopole-based confinement of 
any other fundamental matter. Clearly, the obtained temperature Tgj^^ 
is exponentially larger than (^^'^ just as Tbkt is, that fully validates the 
idea of dimensional reduction. Note also that, at Nf ^> msiK \l, AttT / g'^\ , 
yJlJj) -^ m^ , which means that such a large number of zero modes shrinks 
the molecule to the minimal admissible size. In conclusion, notice that, in 
the presence of W-bosons, the influence of heavy fundamental scalar mat- 
ter to the dynamics of the deconfining phase transition has been studied in 
Ref. [56]. 

8.4. 3D GG model and its supersymmetric generalization at 
finite temperature 

The crucial difference in the finite-temperature behavior of the 3D GG 
model [52] from that of 3D compact QED [51] stems from the presence 
in the former of the heavy charged VF -bosons. Although being practically 
irrelevant at zero temperature due to their heaviness, at finite temperature 
these bosons, due to their thermal excitation, form a plasma whose density 
is compatible to that of monopoles. The density of the VF^-plasma can be 
evaluated as follows (see e.g. Ref. [94]): 



^^ = -'^|^/(^^" 



1-e' 



m-e(p)) 



= 6 

A=o 



d'^p 1 



oo 



(27r)2 e(^<P) 



3m_|. [ dzz ^ 3m|, [ ^^^^-m^pz ^ ^^wT (^^J_\ ^-^wP 
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Here, fl stands for the chemical potential, e(p) = a/p^ + rn'^^^, and the factor 
"6" represents the total number of spin states of W~^- and VF^-bosons. We 
have also denoted z = e{p)/mw and taken into account that, since we 
want to stay in the plasma phase of monopoles and the temperatures under 
consideration should not exceed Tbkt, T <C mw- Therefore, up to small 
corrections, pw = ^Il^e''^^^. 

It has been found in [52], further elaborated in [53, 56], etc., and re- 
viewed in [95] that, when the finite-temperature 3D compact QED is pro- 
moted, upon the incorporation of W-bosons, to the finite-temperature 3D 
GG model, the critical temperature Tbkt and the U(l) universality class of 

2 

the deconfining phase transition are changed to Tc = £^ and the 2D-Ising 
universality class, respectively. In what follows, we will survey the main 
ideas of Ref. [52], proceeding further to the analysis of the supersymmetric 
3D GG model at finite temperature [54]. 

The Lagrangian of 3D compact QED, rewritten in terms of the vortex 
operator V = e~*S'™x/2^ jg 

^3D = \{d^xf - 2Ccos(<7™x) = 4- I^M^I' -C[V^ + {V*f] ■ 

In this form it has explicit magnetic Z2-symmetry [50]. At zero tempera- 
ture, this symmetry is spontaneously broken, since {V {x)V* (0)) — > 1. It 
can be shown [50] that the breakdown of the magnetic symmetry implies 
confinement. Since this symmetry is inherited in the 3D compact QED from 
the initial 3D GG model, the deconfinement phase transition, occurring in 
the latter at a certain temperature (equal to Tc), should be associated with 
the restoration of this very symmetry. Already this fact alone indicates that 
the phase transition in the finite-temperature 3D GG model should be of 
the same kind as that of the (Z2-invariant) 2D Ising model, rather than 
the BKT phase transition of compact QED. It has been shown quantita- 
tively in Ref. [52] that this is indeed the case if W-bosons are additionally 
included in the compact-QED Lagrangian. At finite temperature, this can 
be done by noting that W-bosons are nothing but vortices of the X'^^ld 
and they can be incorporated into the dimensionally-reduced Lagrangian, 
^2D = |(9^x)^ ~ 2,^cos(\AK^x)) by adding to it the term — 2;Ucosx (from 
now on, X will denote the dual-photon field in the dimensionally-reduced 
theory). The field x, dual to the field X) is defined through the relation 
id^x = gVP^fiuduX- The fugacity of W-bosons, fi, is proportional to their 
density pw, therefore /U oc mwTe~"^^^. [Note also that an alternative way 
to introduce VF-bosons into £20 has been proposed in Ref. [96]. Within 
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that approach, the above definition of x through x is abolished. Instead, an 
extra interaction of these, now independent, fields of the type i{dxiX){dx2X) 
appears.] 

Owing to the novel cosine term, even in the absence of the monopole 
plasma, the dual photon never becomes massless [52]. Rather, its mass mu 
increases, and the vacuum correlation length d decreases, with the increase of 
the temperature. Consequently, contrary to what we had in the case of the 
inverse BKT phase transition, where the correlation length was becoming 
infinite at T > TbkT) now it never becomes infinite. This result parallels 
the general expectation that, with the increase of the temperature in any 
local field theory, the degree of disorder becomes higher, and the correlation 
length decreases. Indeed, in the absence of monopoles (i.e. at ^ = 0), 



C2B = l:i^f,c^?-2ficos(g^/^(f>), (8.35) 



where (j) = VTx/g- Therefore, mjj{(3) = 2n[3g^ oc mwg^er"^^^ ■ We see 
that rriD grows with the decrease of /3, i.e. with the increase of T, proving 
the above statement. 

Evaluating further the mean squared separation in the VF"'"I^~-molecule 
in the same way as we did earlier for the MM-molecule [cf. Eqs. (8.1) 

2 

and (8.34)], we have {L"^) oc f (Pini'^~2^ . The convergence of this integral 

2 

at T < TxY = §- means that, in the absence of monopoles, with the increase 
of temperature W-bosons pass from the molecular phase into the plasma 
one at T = Txy- In the same case of absence of dynamical monopoles, a 
static MM-pair becomes linearly confined at T > Txy- Indeed, the MM- 
potential, V, is related to the correlation function of two vortex operators as 
(y(x)y*(y)) = e~l^^(^~y>. In the phase where static monopoles are confined, 
one has 

(y(x)y*(y)) 1''"-^°^ e-^"l^-yl , (8.36) 

and the magnetic Z2-symmetry is restored. Let us clarify the origin of 
Eq. (8.36) and get an idea of the mass parameter which appears in a. To 
this end, notice that we are dealing with the situation dual to that of finite- 
temperature 3D compact QED below Tbkt [51]. In the latter case, external 
static quarks, represented by Polyakov loops, are linearly confined. There- 
fore, just as the correlation function of two Polyakov loops in that theory 
decreases with the distance at the inverse soliton mass, so does the correla- 
tion function of two vortex operators in our case. In fact, it can be shown 
that vortex operators create solitons in the theory (8.35). Solitons in this 
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theory carry a unit topological charge, corresponding to the U(l) symmetry 
generated by the topological current j^ = ^^e^udv(l>- In the dimensionally- 
reduced theory, the 3D vortex operator, V{x) = e"*^™-^*^^''^, is replaced by 
y(x) = e~*^^^(^)/^, or equivalently. 



y(x) = exp I -rr^^ii' / dx^dycj) \ . (8.37) 



'K 

Y 

X 

The integration contour here goes either to infinity or to a point where a 
vortex operator of the opposite magnetic charge, 

y 



V*{y) = exp ^e^, I dx^d,4> , (8.38) 

is defined. One can further demonstrate that the operators (8.37) and (8.38) 
have topological charges 1 and —1 and are nothing but the creation operators 
of a soliton and an antisoliton respectively. This fact stems from the relations 
{QVm = {Vm and {QV*m = - (y*(0)), where Q = e^.j^dx^U^) 
is the topological-charge operator with the contour C encircling the origin 
anticlockwise. These relations stem from others, 

0-^(x)y(o)) c^ -^ {vm , (i^(x)y*(o)) .^ ^ (y*(o)) 

(where "~" means "neglecting correlation functions of topological currents 
higher than the two-point one"), that themselves can be obtained from the 
two-point correlation function of topological currents, 

0.W.(0)> = ||(22g^-V);^. 

Owing to the fact that the operator (8.37) [(8.38)] creates a soliton [an- 
tisoliton], the large-distance behavior of the correlation function (8.36) is 
determined by the lightest state with unit topological charge, i.e. by the 
soliton. One can therefore identify (3(7 on the r.h.s. of Eq. (8.36) with the 
soliton mass m^oh therefore a = Tmsoi- 

We have considered above the two idealistic situations, namely when 
either dynamical monopoles or W-bosons are absent. In reality, when both 
are present and interact with each other, the deconfining phase transition 
temperature is neither Tbkt nor Txy, but is rather Tc, which in fact lies 
between the two. The deconfining phase transition takes place when the 
densities of monopoles and W-bosons become equal, i.e. 2^ = p^. Up to 
inessential preexponential factors, this happens when the monopole action 
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5*0 is equal to mwP, i-e. at T = Tc. At this temperature, the thickness of 
the string connecting a W^- and a VF^ -boson is of the order of its length. 
As has been discussed [cf. e.g. after Eq. (3.17)], the thickness of the string 
is the inverse Debye mass of the dual photon, d oc C~^'^. The length of 
the string is clearly of the order of the average distance between W's, that 
is /i^^'^. Therefore, again up to preexponential factors, the thickness and 
the length of the string are equal at T = Tc. Note that these qualitative 
arguments can be formalized by the RG procedure [52]. In the presence of 
W-bosons, the RG equations possess three fixed points. The first two are the 
zero- and infinite-temperature ones. The third fixed point of the RG flow is 
a nontrivial infrared unstable one, T = T^ S, = IJ-, which corresponds to the 
phase transition. [In this point, however, both fugacities become infinite. It 
has been demonstrated in Ref. [96] that a correcting factor (e + 2)/(2e + 1), 
by which Tc should be multiplied, appears if one demands that the RG flow 
of the fugacities should stop when at least one of them becomes of the order 
of unity.] 

Note that an independent condition of the deconfining phase transi- 
tion stems from the coincidence of the scaling dimensions of the operators 
:cos (vi^x) • ^-iid :cos {g^/JScfj : at the critical temperature. The scaling di- 
mensions are equal ^ and ^^ respectively, so that the monopole cosine term 
is relevant at T < TbkT) whereas the cosine term of the W-bosons is relevant 
at T > TxY- The two scaling dimensions become equal when g'^T = g'^(3, 
that yields the critical temperature |^. As we see, this value is indeed equal 
to Tc, up to the factor -, which is of the order of unity and is generated 

2 

by the loop corrections. At the temperature §-, both scaling dimensions 
are equal to unity, therefore both cosine terms in C2T) are relevant at this 
temperature. In fact, in the whole region of temperatures Txy <T < Tbkt, 
both terms are relevant. 

The supersymmetric version of the 3D GG model [54] contains, in addi- 
tion to the dual photon, its superpartner - an adjoint fermion, which we will 
call the photino. This model possesses a discrete parity symmetry, which 
should lead to the masslessness of photinos. At zero temperature, this is 
however not the case, since the parity is spontaneously broken via a nonva- 
nishing photino condensate. Thus, at flnite temperature, one can anticipate 
two phase transitions - one related to the vanishing of the photino conden- 
sate and the other one due to deconfinement. These two transitions could 
either be distinct and happen at different temperatures, or could coincide. 
In this respect, the model is similar to QCD with adjoint quarks, where a 
similar question can be asked about the (non-)coincidence of deconflnement 
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and restoration of discrete chiral symmetry. 

The Lagrangian of the model contains the bosonic fields of the nonsuper- 
symmetric GG model, that are the massless photon, the heavy W vector 
bosons and the massive Higgs field, as well as their superpartners - photino, 
winos and Higgsino. It has been shown in [97] that, just like in the GG model, 
the monopole effects render the dual photon massive, although the mass in 
this case is parametrically smaller, since it is due to the contribution from 
a two-monopole sector, rather than a single-monopole sector as in the GG 
model. The low-energy sector of the theory is described by the supersym- 
metric sine-Gordon model. Its Euclidean action in the superfield notation is 
[We adopt here the notations of Ref. [98], in particular f (f996 = 1 .] 



s = - d^xcfe 



- <PDaDa^ + C cos {gm^] 



(8.39) 



In this equation, the scalar supermultiplet and supercovariant derivatives 
have the form 



2 ada <j9a 

Here, x again denotes the dual-photon field (real scalar), A is the photino 
field, which is the two-component Majorana spinor (A = A^r^), F is an 
auxiliary scalar field, d = 7i5j, and 7 = r. The monopole fugacity C, has 
dimensionality [mass]^ and is exponentially small. In terms of the disorder 
operator, the action (8.39) in component notations can be rewritten, up to 
an inessential constant, as 

S = fd^x [4- \d^V\'- i A9A - ^ {V^+ y*2) AA - ^^^ {V^+ V*^) 
J \_9in 2 2 2 

where C = ^/4. Besides the magnetic Z2-symmetry, this action has an 
additional discrete parity symmetry inherited from the full supersymmetric 
GG action, 

y(xi, X2, X3) -^ iV{-Xi,X2,X3), A(xi, X2, X3) -^ r^A(-Xi, X2, X3) . 

The photino mass term is odd under the parity transformation. Thus, the 
photino can acquire a mass only if parity is spontaneously broken, that is 
indeed the case. The breaking of parity results in a nonvanishing photino 
condensate (AA) ~ g^C^^W and leads to a nonvanishing photino mass m = 
'^9rni (equal to the Debye mass of the dual photon). One can prove [54] that 
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the equality of the dual-photon and the photino masses is preserved on the 
quantum level as well. 

At finite temperatures, one can integrate photinos out in the 
dimensionally-reduced theory. Again including W-bosons, we arrive at a 
theory similar to the one we had in the nonsupersymmetric case. Its action 
reads 

Sd.-r. = J d^x ]^{d^x? - 2ecos (2/^x) - 2^ cos [g^/^(t) , (8-40) 

where ^ oc /3C^ is a certain positive and exponentially small fugacity. The 
crucial difference from the nonsupersymmetric GG model is due to the fac- 
tor "2" in the term cos(2vXx)) which makes the model (8.40) Z4-invariant. 
Similarly to the nonsupersymmetric case, up to higher-loop corrections, the 
condition for the determination of the deconfinement phase transition tem- 
perature reads 4g'^T = g'^/3, that yields Tc = ^. The universality class 
of the phase transition is therefore Z4 rather than Z2 , as it was in the non- 
supersymmetric case. Further, it has been argued in Ref. [54] that the parity 
order parameter, {V'^'j + (V*'^'), vanishes at T > T^. Therefore, the parity 
restoring phase transition takes place at the same temperature Tc as the 
deconfining transition. While this is an interesting phenomenon, it seems 
to be somewhat non-generic. In particular, in 4D gauge theory with adjoint 
fermions there is no reason to expect the deconfining and chiral-symmetry 
restoring phase transitions to coincide. The physical order parameter for 
deconfinement is the 't Hooft loop V [99] , while for chiral symmetry it is the 
fermionic bilinear form AA. In four dimensions, the two have a very different 
nature. While AA is a local field, F is a string-like object. It is thus difficult 
to imagine these two order parameters combining into a single one as it is the 
case in the 3D theory discussed above. The lattice results indeed indicate 
that, at least in the SU(3)-theory, in four dimensions the two transitions are 
distinct [100]. 

9. Summary 

Below we would like to mention once again some specific issues discussed in 
this review. 

• 3D GG model at T = 0, iV > 2 

— niH < 00 ^ N must be smaller than exp [(certain constant) x 
mw/g'^], otherwise the Higgs vacuum is not stochastic; 

— niH = 00 ^ the string tension of the fiat Wilson loop in the fun- 
damental representation is obtained; it possesses an ambiguity in 
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the numerical factor due to the exponentially large thickness of the 
string; 

— the Kalb-Ramond field, which incorporates both monopoles and 
free photons =^ the numerical factor at all the string coupling con- 
stants (in particular, at the string tension) is fixed for an arbitrar- 
ily shaped surface, in the weak-field (low-density) approximation; a 
generalization of the so-obtained theory of confining strings to the 
adjoint case and to A;-strings: -^^^ ~ 2 at iV > 1; -^^ = ^^)^ 
(Casimir scaling) in the low-density approximation; the leading non- 
diluteness correction is derived, such that, at k r^ N ^ 1, it can 
significantly distort the Casimir scaling; 

— 4D-case with the field-theoretical 0-term in the strong-coupling 
regime ^ the string 9-teim.; fundamental, adjoint, and k-case critical 
values of 9, at which crumpling might disappear; a modification of 
the vacuum structure due to the 9-term.. 

• 3D GG model at T > Tj.r. ~ C^ > effects of W-bosons are neglected 

2 

— N = 2, niH < oo =^ -^— may not be larger than a certain function 
of -^^^ otherwise the Higgs vacuum is not stochastic; 

— Higgs-induced corrections to the BKT RG flow in the leading order 
in mjj ; mn itself evolves very slowly in the vicinity of the BKT 
critical point; A^ > 2 =^ unlike the SU(2)-case, the RG invariance 
holds only modulo the approximation ^ajjCos[(gi — qj)b] ~ "^au, 

ij i 

even at niH -^ cxd; 

— in the presence of Nf dynamical fundamental quark flavors, at 

2 

Nf = 1 and ruq = 0, Tbkt = f^! at Nf > 2 and/or niq ^ 0, 
any fundamental matter (including quarks themselves) is deconfined 
at T > Tdj..; 

• W-bosons are taken into account, supersymmetric generalization =^ 
the deconfining and the discrete-parity restoring phase transitions 
occur at the same temperature, |^; the universality class of the de- 
confining phase transition is the same as in Z4-invariant spin models. 

• confining string 

— the presence of the negative-stiffness term forces the introduction 
of a term into the string action that suppresses the formation of 
spikes and, thus, prevents the crumpling of the world sheet. This 
action, in the large-D limit, has an infrared-stable fixed point at 
zero stiffness, that corresponds to a tensionless smooth string with 
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Hausdorff dimension 2; 

— the effective theory describing the infrared behavior of the confining 
string is a conformal field theory with central charge c = 1; 

— at high temperature, the free energy of the confining string goes 
like -F^(/3) oc — 1//3^ with (3 = 1/T and it agrees in sign, tempera- 
ture behavior, and the reality property with the large- A^ QCD result 
obtained by Polchinski. 
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